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Nazev prace: Reguldrm a chaoticke pohyby v magnetosfere cemych der 
Autor: Ondfej Kopdcek 

Katedra(ustav): Astronomicky ustav Univerzity Karlovy 

Skolici pracoviste: Astronomicky ustav Akademie ved Ceske republiky 

Vedouci doktorske prace: doc. RNDr Vladimir Karas, DrSc, Astronomicky ustav 

Akademie ved Ceske republiky, Praha 

Abstrakt: Pusobem silne gravitace v okoU cemych der muze vest k urychlem hmoty. 
Vteto prdci zkoumdme vlastnosti systemu tvofeneho rotujici cemou dirou ve velko- 
rozmerovem uspofddanem magnetickem poll Nabite cdstice v bUzkosti horizontu 
jsou krome silne gravitace ovlivnovdny magnetickym polem a indukovanym polem 
elektrickym. Oprotijiz drive v literature diskutovanym situacim pfiddvdme nekterd 
vyznamnd zobecneni. Magneticke pole v nasem pnpade nemusi byt koaxidlni s ro- 
tacni osou cerne diry, takze system ztrdci osovou symetrii. Krome toho pfedpo- 
kldddme translacm pohyb cerne diry s obecnym smerem i rychlosti. Ukdzeme, ze 
dfky tomu dochdzi k novym efektum. V komplikovane struktufe vysledneho mag- 
netickeho pole pozorujeme v ergosfefe jeho rychle prostorove zmeny provdzene 
vznikem nulovych bodu, ktere dokazuji, ze gravitacm pusobem rotujiciho zdroje 
muze podnecovat rekonekci magnetickych silocar Ddle zkoumdme dynamicke vlast- 
nosti nabitych cdstic vystavenych pusobem tohoto typu poll. Pfedevsim se zajimdme 
o pfechody mezi reguldmim rezimem a deterministickym chaosem, ke kterym do- 
chdzi V zdvislosti na volbe parametru. Pfi numerickem zkoumdni cdsticove dy- 
namiky aplikujeme v kontextu obecne relativity zatim nepouzitou metodu rekurencni 
analyzy. 

Klicova slova: obecnd relativita, kompaktni telesa, astrofyzikdlni korona, deter- 
ministicky chaos 
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Title: Transition from regular to chaotic motion in black hole magnetospheres 
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Abstract: Cosmic black holes can act as agents of particle acceleration. We study 
properties of a system consisting of a rotating black hole immersed in a large- 
scale organized magnetic field. Electrically charged particles in the immediate 
neighborhood of the horizon are influenced by strong gravity acting together with 
magnetic and induced electric components. We relax several constraints which 
were often imposed in previous works: the magnetic field does not have to share 
a common symmetry axis with the spin of the black hole but they can be inclined 
with respect to each other, thus violating the axial symmetry. Also, the black hole 
does not have to remain at rest but it can instead perform fast translational motion 
together with rotation. We demonstrate that the generalization brings new effects. 
Starting from uniform electro-vacuum fields in the curved spacetime, we find sep- 
aratrices and identify magnetic neutral points forming in certain circumstances. 
We suggest that these structures can represent signatures of magnetic reconnection 
triggered by frame-dragging effects in the ergosphere. We further investigate the 
motion of charged particles in these black hole magnetospheres. We concentrate 
on the transition from the regular motion to chaos, and in this context we explore 
the characteristics of chaos in relativity. For the first time, we apply recurrence 
plots as a suitable technique to quantify the degree of chaoticness near a black 
hole. 

Keywords: general relativity, compact objects, astrophysical coronae, determinis- 
tic chaos 
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Chapter 1 

Introduction 



1.1 Astrophysical black holes 

Nowadays it is a consensus that black holes are a vivid part of physical reality. 
Astrophysicist's attitude toward the possibility of a real existence of completely 
gravitationally collapsed bodies, i.e. black holes (BHs), has undergone profound 
changes since 1 935 w hen sir Eddington commented recent theoretical results of 



Chandrasekharl (119351 ). suggesting that black hole could be the endpoint of star 
evolution, by these words: "I think there should be a law of Nature to prevent a star 
from behaving in this absurd way". At present, due to discoveries and observations 
made in recent decades, it seems highly probable that there is actually no such 
power which would save a heavy star from collapsing into the black hole. 

Until the late 1960's there had not been much progress made in this field. Al- 
though Karl Schwarzschild had given an exact solution to Einstein's equations for a 
spherically symmetrical source (Swcharzschild black hole in the case of collapsed 
body) already in 1916 (just one year after publishing the theory of relativit y), it 
was n ot believed that such objects really exist. Almost a half century later, i Kenj 



(119631 ) gave a new exact solution to Einstein's field equations describing the geom- 



etry around a rotating compact object, which naturally appears to be more relevant 
in astrophysical context than prior non-rotating Schwarzschild solution. But this 
progress still would not become a concern of astrophysicists unless there had not 
been made important observational discoveries during 1960's. 

Astronomers discovered sources of radio waves with highly redshifted (with 



redshift factor z ^ 0.1 — 6) spectrum ( Schmidt . 19631 ) which proved them (accord 



ing to Hubble's law) to be very distant from us 0.3 — 3 Gpc ). Combined with 
observed visual magnitude of these objects - called quasars (quasi- stellar radio 
sources) - there comes a conclusion of their huge power output of about 100 times 
that of the total luminosity of average galaxy (Lquasar ^ 10^^ — 10"^^ W). Quasars 
also appeared to emit significantly in X-ray and even in gamma part of the spec- 
trum. 

Quasars were not the only peculiar objects observed during those decades. Ac- 
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cording to the type of their spectrum and luminosity time dependence of those 
"new" objects they were classified as Seyfert galaxies (being observed and stud- 
ied since the 1940's actually), blazars, or radio galaxies (RGs) with all of them 
consisting of further subclasses. They are all characterized by extraordinarily high 
luminosity coming from a small volume 10~^ pc^). As there is no consistent 
way to explain the mechanism of those energy sources conventionally (i.e. regard- 
ing stars and therein running nuclear synthesis as the most efficient energy source 
in the universe), attention was turned to hypothesis employing strong gravitational 
fields considering compact objects and subsequently the black holes. Now it is gen- 
erally believed that above mentioned phenomenons are of the same origin which 
was given the name active galactic nuclei (AGN). Galaxies with active galactic 
nuclei and subsequently also those with non-active nuclei are suspected of hosting 
a supermassive black hole (MBH) of M ^ 10^ — l(f Mq in their centres, e.g. M87 
- AGN of ^ 10^ or Sagittarius A* - nonactive nucleus of our Galaxy - with 
M^4Ax 10^ Mq (lGenzel etal.Ll20ld ). 





Figure 1.1: Same physical mechanism is supposed to operate in different types of 
cosmic objects. Central black hole is accreting matter which forms an accretion 
disk. Collimated jets ar e launched from the central part of the system. Illustration 
credit: iMirabell (l2007bb . 



The main reason for such an assumption is that it provides a clarification of 
the observed luminosity of active galaxies. When the Kerr black hole is employed 
to model the properties of AGN we conclude that the matter accreted from its 
vicinity would form an accretion disk which becomes heated by accretion process 
and subsequently emits radiation of various types (depending on temperature and 
many other p r operti es of the disk). Doing so there could be as much as ^ 40% 
(iMisner et al. I 09731) p. 885) of accreted material rest mass turned into emitted 
radiation which is considerably more than in the case of thermonuclear synthesis 
of helium (^ 7% of the rest mass). 
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Whether we regard observed AGN as quasar, radio galaxy, blazar or Seyfert 
galaxy depends mostly on the spatial orientation of the surrounding galaxy which 
acts as the shie ld for some parts of the spectrum while being a source in some 
other parts (see lUrry & Padovani Il995l . for a review of unified schemes of radio 
sources). 

Tolman-Op penheimer-Volkoff (TOV) limit sets the upper bound to the mass 



of neutron star (lOppenheimer & Volkofn . Il939h. Today's estim ates of the value of 



TOV limit range among ^ 1.44 -3Mq dLattimer & Prakashl . 12004) and the un- 
certainty is due to unknown equation of state. Pressure of degenerate neutron gas 
in the neutron star reaching this limit cannot oppose gravitational pressure and in- 
evitably collapses into the black hole (while collapse is being accompanied with 
gamma ray burst, GRB). Stellar mass black hole (M ^ 4 — 15 Mq) should thus rep- 
resent the final state of stellar evolution of every heavier star. To be observed it 
needs to be a part of the binary system. Under certain circumstances black hole 
accretes matter from its companion star (flowing via Lagrange point when bound- 
aries of a Roche lobe are exceeded). Similarly like in the supermassive black hole 
case there establishes accretion disk emitting in X-ray part of the spectrum. Due to 
recent observation from the satellites Chandra and XMM-Newton equipped with 
high resolution X-ray detectors we now have many stellar mass black hole can- 
didates in our Galaxy (e.g. Cygnus X-1). The key to distinguish accretion disk 
of a compact object (typically neutron star) from the one surrounding the black 
hole is the nonexistence of the surface in the black hole case. Thus irregular flares 
of gamma rays accompanying thermonuclear reaction of accreted material on the 
compact object surface are not detected. 

Intermediate-mass black holes (IMBHs) are those with M ^ 10^ — W"^ Mq. 
Existence of IMBHs is still uncertain although a number of candida tes was iden- 
tified . IMBH might be possibly formed in globular star clusters ( Maccarone et al.L 



2007|). Ultra-luminous X-ray sources (ULXs) in close galaxies are suspected to be 



powered by IMBH. An U LX source HLXl l ocated on the edge of spiral galaxy 
ES0243-49 was claimed (Farr ell et all . l2009l ) to host a n IMBH of oy er 500 M®, 
though the interpretation was recently disputed by Soria et al.l (1201 ih . A stellar 
complex IRS 13E residing close to our galactic center Sagittarius A* is yet another 
candidate for the IMBH system dSchoedel et al.l . l2005b . 

Primordial black holes (PBHs) are hypothetical objects whose origin differs 
fundamentally from all above mentioned types. They have not been established 
by gravitational collapse of stars or other astronomical bodies since they might 
have been born from density fluctuation during early stages o f the evolution of 
the U niverse. On the theoretical grounds it has been argued / Carr & Hawking . 
19741 ) that PBHs of mass from 10~^ g upwards migh t exist in th e present Universe. 



Updated constraints on PBHs were given recently by lCarr et al. 1(2010). Primordial 
black holes could be detected due to the Hawkin g radiation which black holes emit 
according to the quantum relativity ( Page . 19761 ). 

Hawking radiation has ordinary thermal spectra and energy it radiates away 
goes on account of the mass - the black hole evaporates. Predicted evaporation 
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times for stellar mass black holes exceed present age of the universe by many 
orders but for primordial black holes with the mass M ^ 10^^ kg evaporation time 
approaches its current age. Rate of evaporation escalates as the mass decreases 
(power output P tjt). Endpoint of the primordial black hole existence should 
thus be explosive Mawkingl . 1 19741 ) and is believed to be accompanied with loud 
GRB. Search for these GRB signatures is one of the key scientific objectives of 
current mission of Fermi Gamma-ray Space Telescope which is operating since 
2008. 
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Figure 1.2: Evidence for fast motion of the emerging jets in the microquasar GRS 
1915+105 (observed at radio frequency of 8.6 GHz) and in the quasar 3C 279 (at 
22 GHz). Synchrotron emission has been reported in infrared wavelengths and, in 
some cases, even up to X-ra ys, implying the presence of electrons in the jets with 
TeV energies. Figure credit: iMirabel (l2007al ). 



1.2 Magnetic field: trigger for accretion and outflow 

Recent observations of microquasars, pulsars, gamma-ray bursts indicate that the 
astrophysical jets play an important role everywhere (not only in the case of AGNs). 
There is plenty of observational evidence suggesting that the initial acceleration of 
jets takes place very near black holes (or other compact object) and proceeds via 
electromagnetic forces. Jets and accretion disks in the vicinit y of compact objects 



probably create symbiotic magnetically driven system (e.g. iFalcke & Biermann 



19951). 

The current promising model of the dynamics (i.e. launching, accelerating 
and collimating) of the astrophysical jets is based on the magnetohydrodynamics 
(MH D). The results of the sim ulations employing general relativistic MHD equa 



tions (IKrolik & Hawleyl . l20ld ) correlate with observations of M87 (IJunor et al.l . 



19991 ) where the formation and the coUimation of the jet were analyzed. 
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Moreover, the 3D relativistic MHD simulations carried out by lHawley & Krolik 



(120061 ) demonstrate clearly the essential role which the accretion disk's coronae 
play in the collimation and acceleration of the jet. Indeed the dominant force accel- 
erating the matter outward in a given numerical model originates from the coronal 
pressure. Regions above and below the equatorial plane become dominated by the 
magnetic pressure and large-scale magnetic fields may also develop by the dynamo 
action. 

Recent numerical relativistic study by iRezzoUa et al.l (120111 ) reveals the for- 
mation of ordered jet-like structure of ultrastrong magnetic field in the merger of 
binary neutron stars. Such system thus might serve as an astrophysical engine for 
observed short gamma ray bursts. 

In accretion models the magnetic field was also employed - the so called mag- 
netorotational instability (MRI) must operate in accr etion disc, gen erating the ef- 
fective viscosity necessary for the accretion process (lBalbusl . ll99ll ). Magnetic re- 
connection is likely to be responsible for rapid flares, which are observed in X- 
rays. Finally, Fara day rotation measureni ents suggest that tangled magnetic fields 
are present in jets (IBegelman et al.l . ll984l ). 

Observations of the Galactic Center (GC) reveal the presence of another re- 
markable large-scale magnetic structure - nonthermal filaments (NTFs). NTFs 
cross the Galactic plane and their length reaches tens of parsecs while they are 
only tenths of parsec wide. The strength of the magnetic field withi n the NTF may 
appro ach ^ 1 mG while the typical interstellar value is ^ lOjuG (ILaRosa et al.L 



20041) . Initially, it was thought t hat NTF s trace the pervasive poloidal magnetic 



field present throughout the GC ( Morrisl . [l990). Later, however, it became ap- 
parent that the str ucture of t he ma gnetic field in the central region of the Galaxy 
is more complex dFerriere . 2010 ). See fig. 11.31 for the snapshot of the GC at 
90 cm (330 MHz) which shows the NTFs clearly. 

Overall it is quite likely that electromagnetic mechanisms play a major role 
and operate both near supermassive black holes in quasars as well as stellar-mass 
black holes and neutron stars in accreting binary systems (see fi gure s [TT] and 11.21) . 
Besides that a faint magnetic field is present throughout the interstellar medium, 
being locally intensified in NTFs. 



Electro-vacuum fields 

The survey of the vacuum electromagnetic (EM) fields may be regarded as the fun- 
damental starting point in studying the dynamics of astrophysical plasma. If the 
examination of the fields itselves represents the initial step in a given direction, at 
succeeding stage we would consider the motion of a non-interacting test particles 
exposed to these fields. In this work we will be dealing with both issues. Structure 
of a particular astrophysically motivated EM field emerging in the vicinity of rotat- 
ing black hole will be studied in detail. Subsequently we shall discuss the motion 
of charged particles exposed to the field representing a special case of a general 
solution explored before. Primarily we concern ourselves with the stable orbits oc- 
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Figure 1.3: Milky way's Galactic Center penetrated by the narrow nonthermal fil- 
aments (NTFs). The strength of ordered magnetic field may approach ^ 1 mG 
within NTFs. Length of NTFs reaches tens of parsecs. Inner region of 0.8° x 1.0° 
is shown at waveleng th 90cm (330 MH z ). Snapshot was take n by the Very Large 
Array (VLA). Credit: iNord et all (l2004l) : lLaRosa et all t004 ). 



cupying off-equatorial potential lobes. Particles on these orbits are relevant for the 
description of astrophysical corona comprising of diluted plasma residing outside 
the equatorial plane in the inner parts of accreting black hole systems. 

Gaseous corona is supposed to play a key role in the formation o f observed X - 
ray spectra of both active galactic nuclei (AGNs) and microquasars (lDonel . l200lh . 
Power law component of the spectra is believed to result from the inverse Compton 
scattering of the thermal photons emitted in the inner parts of the disk. Relativistic 
electrons residing in the corona serve as a scatterers in this process. Their dynamic 
properties (e.g. resonances) thus shall have imprint on the observed spectra. 

The role of magnetic fields ne ar strongly gravitating objects has been subject of 
many investigations (e.g. IPunsM 12008. ). They are relevant for accretion disks that 
may be embedded in large-scale mag netic fields, for example when the accretion 
flow penetrates close to a neutron star ( Lipunov et al. . 1992 : Mestei 1999b . Outside 
the main body of the accretion disk, i.e. above and below the equatorial plane, the 
accreted material forms a highly diluted environment, a corona, where the density 
of matter is low and the mean free path of particles is large in comparison with the 
characteristic length-scale, i.e. the gravitational radius of the central body, Rg = 
GM/c^ ^ \ .5{M /Mq) km, where M is the central mass. The origin of the coronal 
flows and the relevant processes governing their structure are still unclear. In this 
context we discuss motion of electrically charged particles outside the equatorial 
plane. 
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Regular and chaotic dynamics 

Accretion onto black holes and compact stars brings material in a zone of strong 
gravitational and electromagnetic fields. We study dynamical properties of motion 
of electrically charged particles forming a highly diluted medium (a corona) in the 
regime of strong gravity and large-scale (ordered) magnetic field. 

We start our discussion from a system that allows regular motion, then we 
focus on the onset of chaos. To this end, we investigate the case of a rotating black 
hole immersed in a weak, asymptotically uniform magnetic field. We also consider 
a magnetic star, approximated by the Schwarzschild metric and a test magnetic 
field of a rotating dipole. These are two model examples of systems permitting 
energetically bound, off-equatorial motion of matter confined to the halo lobes 
that encircle the central body. Our approach allows us to address the question of 
whether the spin parameter of the black hole plays any major role in determining 
the degree of the chaoticness. 

The both dynamic systems may be regarded as different instances of the orig- 
inally integrable systems which were perturbed by the electromagnetic test field. 
Complete integrabilit y of geodesic motion of a free particle in Schwarzschild space- 



time is easy to verify (IMisner et al.Lll973h . To some surprise it was later found that 



also free particle motion in Kerr spaceti me and even the charged particle motion in 



Kerr-Newman is completely integrable ( Carted . 1 19681 ) since separation of the equa 



tions of motion is possible as there exists additional integral of motion - Carter's 
constant L . Trajectories found in such a systems are merely regular. 

In the non-integrable system, however, both regular and chaotic trajectories 
may coexist in the phase space. Standard method of a qualitative survey of the 
non-linear dynamics is based on the construction of Poincare surfaces of section 
which allow to visually discriminate between the chaotic and regular regime of 
motion. 

On the other hand quantifying the chaos by Lyapunov characteristic exponents 
(LCEs), as its standard and commonly used indicator, becomes problematic in the 
general relativity (GR) since LCEs are not invariant under the coordinate transfor- 
mations. Besides that the usual method of computing LCEs involves evaluation of 
the distances between the neighboring trajectories which becomes intricate in GR. 
Alth ough there are opera tional workabouts to partially overcome these difficulties 



(e.g. IWu & Huangl . |2003|) the need for a consist ent treatment is apparent. Perhaps 



the geometrical approach suggested recently by lStachowiak & Szydlowskil (120 1( 
could eventually provide a covariant method of the evaluation of the Lyapunov 
spectra in GR. 

In this context w e adopt a different to ol to investigate the dynamic system - 



Recurrence Analysis (IMarwan et al.L l2007h . To characterize the motion, we con- 



struct the Recurrence Plots (RPs) and we compare them with Poincare surfaces of 
section. We describe the Recurrence Plots in terms of the Recurrence Quantifica- 
tion Analysis (RQA), which allows us to identify the transition between different 
dynamical regimes. We demonstrate that this new technique is able to detect the 
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chaos onset very efficiently, and to provide its quantitative measure. The chaos 
typically occurs when the conserved energy is raised to a sufficiently high level 
that allows the particles to traverse the equatorial plane. We find that the role of the 
black-hole spin in setting the chaos is more complicated than initially thought. 

1.3 Structure of Thesis 

The thesis is organized as follows. Most of its contents are contained in chap.O We 
begin with technical preliminaries, namely in sec. 12.21 we introduce several alterna- 
tive definitions of electric and magnetic vector fields. In sec. 12.31 we give explicitly 
the components of the electromagnetic tensor F^y describing the field around the 
Kerr black hole drifting in the arbitrary direction through the asymptotically ho- 
mogeneous magnetic field which is generally inclined with respect to the rotation 
axis of the BH. Choice of the observer's frame is discussed in sec. 12.41 Structure 
of the both electric and magnetic fields is explored in detail in sec. 12.51 First in 
sec. 12.5.1] we revisit the issue of the magnetic expulsion (Meissner effect) which is 
present in the case of aligned field. Subsequently we study field structures in the 
case of inclined field. Finally, in sec. 12.5.2] we introduce a drift of the black hole 
and explore the both electric and magnetic fields emerging in this general setup. 

In the remaining part of the chap. [2] we shall deal with the dynamics of the 
charged test particles exposed to the test fields analyzed in the previous sections. In 
sec. I2.7l we review the equations of particle motion, which we then integrate to ob- 
tain trajectories. In sec. l2.8l we introduce the basic properties of Recurrence Plots. 
Sec. 12.91 analyses the motion around a Kerr black hole endowed with a uniform 
magnetic test field. We employ Poincare surfaces of section and Recurrence Plots. 
The two approaches allow us to show the onset of chaos in different, complemen- 
tary ways. We examine the motion in off-equatorial lobes, pay special attention 
to the spin dependence of the stability of motion, and we notice the emergence 
of 'potential valleys' that allow the particles to escape from the equatorial plane 
along a narrow collimated corridor. Analysis of the off-equatorial motion around a 
magnetic star is presented in sec. 12.101 We consider a dipole-type magnetic field, 
which sets different limits on the off-equatorial range of allowed motion of charged 
particles. It also defines different regimes of chaoticness in the comparison with 
the uniform magnetic field. 

Finally, results of the analysis are summarized in chap. [51 Chapter |4]represents 
a brief outlook to the future as it specifies several topics suggesting the direction in 
which our research could continue. 

In Appendix we list conversion factors between SI and geometrized units. 
Scaling of the quantities by the central mass M is also discussed. Appendix |B]gives 
detailed comparison of numerical integrators applicable to our system. Benefits 
of symplectic routines are discussed therein. In Appendix [C] we introduce our 
program tool vfexplorer2 which provides a simple graphic user interface (GUI) for 
the effective exploration of a complex vector fields. 
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M ost of the scien t ific con tents presented in the Thesis were published in the 
papers iKopacek et alJ (I2010al") and iKaras & Kopacek (l2009h . Investigation of the 
particle motion in the astrophysical corona was pr eceded by the study of the topol- 
ogy of off-equatorial potential lobes performed by iKovaf et all (|20l d). Discussion 
of the charged particle motion in such lobes was also a subject of the contribution 
Kopacek et all (l2010bl ). Initial steps of the investi gation of electro magnetic fields 



around drifting Kerr black hole were described by IKopacek 

At different phases of my doctoral studies I took part in several research projects. 
Namely I acknowledge support from the doctoral student program of the Czech 
Science Foundation (project No. 205/09/H033), Plan for European Cooperating 
States of the European Space Agency (ESA PECS 98040) and two projects of the 
Grant Agency of Charles University (GAUK 119210/2010 and SVV-263301). 



Chapter 2 



Regular and chaotic motion 
in black hole magnetospheres 



2.1 Electromagnetic field around drifting Kerr black hole 



In this section we construct a test field solution describing the electromagnetic field 
around a Kerr black hole which is drifting in an arbitrary direction with respect to 
the asymptotically uniform magnetic field with the general orientation with respect 
to the rotation axis. Definition of the electric and magnetic field intensities is dis- 
cussed as well as the choice of the physical observer. Structure of the field is 
explored in detail. 

Electromagnetic (EM) test field solutions of Maxwell equations in curved space- 
time play an important role in astrophysics since we can usually suppose that as- 
trophysically relevant EM fields are weak enough, so that their influence upon 
background geometry may be neglected. 

We are interested in the solutions describing an originally uniform magnetic 
field under the influence of the Kerr black hole. Since the Kerr metric is asymptot- 
ically flat, this EM field reduces to the original homogeneous magnetic fiel d in the 
asymptotic region. First such a test field solution was given by lWaldl (Il974r) for the 
special case of perfect alignment of the asymptotically uniform magnetic field with 
the symmetry axis. Using a different approach of Newman Penrose formalism a 
more gen eral solution for an arb itrary orientation of the asymptotic field was in- 
ferred by iBicak & Dvof^ (Il980b . We use their solution to construct the EM field 
around the Kerr black hole which is drifting through the asymptotically uniform 
magnetic field. 

Such generalized setup shall represent an astrophysically relevant model. In 
the actual accreting BH system the misaligned ordered field may arise if the accre- 
tion disk is inclined with respect to the rotation axis of the BH and the Bardeen- 
Petterson effect does not operate to align the axes. Such misaligned a ccretion was 
observed in some numerical simulations (e.g. [Rockefeller et al.l . l2005b . A possible 
scenario how the BH could receive substantial velocity with respect to its accre- 
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tion disk is suggested by the simulations of a merger process (e. g. iRezzoUal . l2009l : 



Gonzales et aLLl2007h . 



In the context of black hole mergers iLyutikovl (1201 Ih recently studied electro- 
dynamic properties of the simplified model consisting of the Schwarzschild black 
hole in the uniform transversal motion with respect to the homogeneous magnetic 
field. The author considers the interaction of the field with the plasma generated 
from the pair production due to the vacuum breakdown. Analyzing the resulting 
situation in the force-free approximation, Lyutikov concludes that observed elec- 
trodynamic properties resemble in many aspects the pulsar magnetosphere. 



2.2 Lines of force 

Kerr me tric in Boyer-Linda uist coordinates = r, 9, (p) may be expressed as 
follows dMisner et all 1 19731) : 

(i/ = - ^ [dt -asinQdifif + [(r^ + a^)<icp - adtf + ^ dr^ + LdQ^, (2.1) 



where 



r^-2Mr + a^ Z = + a^cos^9. (2.2) 



We stress that geometrized units G = c = ^ = = 1 are used throughout the text. 
See Appendix lAl for details. 

For the sake of reference it might be useful to list both covariant and contravari- 
ant components of the Kerr metric explicitly 

a^sin^9-A Z Asin^9 

gtt = , <?^^^^' <?ee = Z, g(p(p = — - — , (2.3) 

a sin^ 9 . ^ 9 2\ 
8t(^ g(^t — — (A - r -a j , 

^ AZ' ^ Z' ^ Z' ^ Z Vsin29 A^ ' ^ ^ ^ 

denoting A = (r^ + a^Y ~ ^ sin9^A. Following relation proved useful in the cal- 
culations: AZ = A - lMr[r^ + c^). 

Radial metric component grr and contravariant components g^^ and 
are singular at A = which defines the outer (+) and the inner (-) horizon r± — 
M ± VM'^ — a^ of the black hole. By querying the curvature scalars one finds 
that the singularities present at r± are purely coordinate singularities rather than 
physical ones. Both horizons merge at r = M in the case of extreme Kerr black 
hole a = M. For a > M the horizon disappears and the central singularity becomes 
naked. In the following, however, we will be concerned with the cases a<M only. 
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Dragging of the inertial frames caused by the rotation of the source may be 
characterized by the quantity Q. = —gt(p/g^>^> = 2Mra/A. Coordinate angular ve- 
locity d(p/dr = u^/u^ ("angular velocity as measured at infinity") of the observer 
freely falling from the rest at infinity reads (see eq. 12.531) . 

In the classical electrodynamics we spontaneously define electric and magnetic 
lines of force as a field lines of three-vectors E = {Ej^^Ey^E^ and B = {B^^^By^B^ 
which form a solution of Maxwell equations for a given problem. Definition con- 
dition of the field line r = {x{s)^y{s)^z{s)) parametrized by the parameter s is such 
that its tangent vector ^ is at each point parallel to the vector field itself. Thus for 
(e.g. electric) lines of force we obtain equation 

dr ^ dr dv dz 

— x^ = 0^ — = ^ = — , (2.5) 
dtS" E^j^ /^'^ 

where x stands for the ordinary cross product of two vectors. 

If we recall Lorentz relation specifying the force felt by a particle with electric 
test charge qe and hypothetical magnetic monopole test charge qm moving with the 
velocity v in the external fields E, B 

F = qe{E + vxB)+ qyn{B - v x £), (2.6) 

we conclude that in a given reference frame we can identify the electric intensity 
E with the force felt by the unit electric charge and magnetic induction B with the 
force felt by the unit magnetic monopole charge provided that these charges are 
static in a given reference frame. 

In the covariant language of general relativity the Lorentz force felt by the test 
particle of mass m and the electric charge qe or the magnetic monopole charge q^ 
is expressed as follows: 

dT dx m 

Dw^ Au^ u ^ Q ^m*^u V /o o\ 

= — ;— = ——^r'^nu'^u^ = E'^yU^ (2.8) 

dT dx m 

where i stands for the proper time, a is particle's four-acceleration, u its four- 
velocity and are Christoffel symbols. F^v is electromagnetic tensor and *F^v 
its dual. 

In the close analogy with the classical case it appears natur al to define coor 



dinate components of magnetic and electric fields as follows faanni & RuffiniL 



19731) 



= E^^E^^u\ (2.9) 

However, the demand on the test charge to be static in a given reference frame, i.e. 
to have 4-velocity of a form = (w^ 0,0,0), becomes highly problematic when 
dealing with the non-static geometry of the Kerr spacetime. It is well known that 
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in this case no physical observer may remain static inside the ergosphere whose 
boundary is defined by 

rs = M + VM^-a^cosQ^. (2. 10) 

Since the region near the horizon is typically of a big interest when investigating the 
effects of a strong gravity we do not content us with being limited to the region out- 
side the ergosphere. To this end we generalize the above definition eq. (12.91) in such 
a way that we allow general four-velocity of the test charge = {u^ ,u^, u^^u^). In 
the classical analogy this would mean not insisting on v = when defining the 
vector fields E and B using the classical form of the Lorentz force equation (12.61) . 
Nevertheless all the cases we shall discuss actually have = therefore we may 
write explicitely 

B^ = -g^^ {%,u' -^%^u^) , (2.11) 
^ _^ee ^y^^^t +*FecpW^) , (2.12) 

B^^- [^W ^y^^^t ^^f^y^ ^g^t (*77^^^^+*/7^^^9)] . (2.13) 

Equation for the lines of force eq. (12.51) takes following form in the poloidal plane: 
dr/d9 = B^/B^ and similarly in the equatorial plane: dr/d(p = B^/B^. To obtain 
analogical expressions for the components of electric field E^, E^, E^ we just need 
to omit the minus sign at the beginning of each term and replace *F^v components 
with F^v in eqs. ([rTT1) - (l2l3l) . 

Another issue connected with the definition (12.9b is a question of the normal- 
ization of the basis vectors. We notice that components of magnetic and electric 
fields B^ and E^ are provided in the canonical basis of Boyer-Lindquist coordinates, 
namely in the basis of vectors ^ which are not normalized. This prob- 

lem may be seemingly easily overc ome by expressing "physical components" of 
the fields^ (iHanni & Ruffinil.ll973b as follows 



^physical = Sign(50 = (2.14) 

^physical = sign(5^) ^^^b\ (2. 15) 

physical = sign(^^) = sign(^^) ^ {%^u^-^F^x), (2. 16) 

and analogically for physical components of electric field ^^^fyjicai- 

However, in the case of Boyer-Lindquist coordinate system which is singular 
at the horizon (namely the radial metric component grr diverges here) using above 
defined physical components turns out to be a controversial choice since it intro- 
duces this divergence directly into the radial component of electric and magnetic 
fields. 



^The adjective physical is sometimes used to refer to the quantities measured in the frame attached 
to a given physical observer. We stress that our usage of the term differs as our physical components 
are expressed in the coordinate frame. 
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To fix this problem we suggest to define "renormalized components" of the 
fields as follows 

^renormalized ~ ^ ^renormalized ~ ^ (2.17) 

Prenormalized"^ ^0' ^renormalized ~ ^ ^0' (2.18) 
Prenormalized = ^(p 5 ^renormalized =^^^95 (2.19) 

where /z^ = 1, /ze = r, /Z(p = rsin9 are ordinary Lame coefficients of spherical co- 
ordinates in flat space. Correspondence between physical and renormalized com- 
ponents lies in the asymptotic region where they become identical since the Kerr 
spacetime is asymptotically flat and Boyer-Lindquist coordinates r,9,(p asymptot- 
ically turn into spherical coordinates describing the spatial part of Minkowskian 
spacetime. Explicitly given the asymptotic behaviour of the metric coefficients is 

^ 1 = hj, gee = hl and gtptp sin^ 6 = and gt^ 0. 

Going even further in "operational flattening" of the curved background one 
could eventually identify the components of magnetic and electric fields directly 
with the components of the electromagnetic tensor F^y. Such an identification 
is fully justified in any local Lorentz frame where F^v takes the form given by 
eq. (12.251) . Nevertheless although the asymptotics of Kerr background is flat, its 
Boyer-Lindquist coordinate basis 1 -forms are not normalized and normalization 
factors hr, Hq and /Z(p must be used accordingly. We define asymptotically moti- 
vated (AMO) components of the fields 

^AMO J^^' ^AMO (2.20) 



^AMO — T~r^ ^AMO — (2.21) 



o(p _ fie — n ii\ 

^^^^ " hM ^''^ " K ' 



AMO components asymptotically coincide with the physical and renormalized 
components provided that the test charge is static. Although AMO components 
do not allow for a direct physical interpretation when applied outside the asymp- 
totic region, it may still be useful to explore them since they do not employ any 
particular observer (four- velocity of a test charge) in their definition and therefore 
it may be easier to acquire intuitive insight into the nature of a given EM field. 

Nevertheless a consistent way to define the electric and magnetic fields should 
provide obvious physical interpretation of the observables measured by a certain 
physical observer at any distance from the center. We let such an observer with 
four- velocity equipped with the orthonormal tetrad e^^^^ measure the Lorentz 
force of eq. (12.91) using his tetrad basis. Tetrad components of the vector fields 
determining desired lines of force are given as the spatial part of the projection 

B« = Bu) = -e^i *FW = -etlXyel, = -*Fu)u) , (2.23) 



= £(,■) = e^^FW = = /^(OW. (2-24) 
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where e^^^ are 1 -forms dual to the tetrad vectors e^^^y Lowering/raising of the 
spatial tetrad indices does not matter since the tetrad is supposed to be orthonormal: 

^(,<)(v) ='n(^)(v). 

For the sake of the future reference we review both tetrad components of EM 
tensor ^(a)(p) and its dual '^(a)(|3)- We remind that such an interpretation of these 
tensors is possible only in any local Lorentz frame which may be attached to any 
physical observer, not necessarily inertial. Magnetic and electric test charges which 
are used to measure the fields are at rest in this frame. 



'(«)(P) 



/ 



B 



He) 

(9) 



-B, 



(9) 



-B 

B< 



(e) 



(2.25) 



(a)(P) 




-B 



Bi 



(e) 



He) 



^(e) 


-£(r) 



^(9) A 
/ 



(2.26) 



Several ways to define electric and magnetic field have been described in this 
section. Coordinate components of the Lorentz force felt by the unit electric/mag- 
netic test charge (eq. 12.91 ) provide a natural generalization of the classical defini- 
tion. In the case of the Kerr background, however, we cannot insist upon the usage 
of the static test charges in the definition as no static observers may penetrate in- 
side the ergosphere. Since the structure of the fields in the vicinity of the horizon 
is usually of the utmost interest we have assumed more general four-velocity of 
the test charge in the form m'' = {u\u'',0,u'^). The Lorentz force may be projected 
onto the tetrad basis attached to the test charge. 

Nevertheless if we remain in the Boyer-Lindquist coordinate basis we note that 
the fields do not come with the proper dimension because the coordinate basis vec- 
tors are not normalized. To correct this in a rigorous manner we define physical 
components of the fields by eqs. (I2.14l) - (l2.16b . It appears, however, that this def- 
inition amplifies the effect of the coordinate singularity at the horizon. Thus we 
suggested to define renormalized components (eqs. l2.17H2.19t as an useful ap- 
proximation which is less problematic near the horizon. Far more approximative 
(in the sense of treating the curved background as flat) are asymptotically motivated 
AMO components (eqs. l2.20H2.22b which in fact map onto the flat surface di- 
rectly. Brief discussion regardin g the definition of the lin es of force of electric and 
magnetic fields was also held by IB leak & Dvorak (Il980b . 



2.3 Electromagnetic field 

Stationary and axisymmetric test field solutions to the Maxwell equations on the 
Kerr background traditionally attract attention by both relativistic theoreticians and 



CHAPTER 2. MOTION IN BLACK HOLE MAGNETOSPHERES 



16 



astrophysicists. The latter are usually concerned with the solutions describing some 
potentially realistic electrody namic scenario. I n particular w e mention test fields of 
axisy mmetric cur r ent loops dPetterson ^ 1975 : Moss . 2011 ) and uniform magnetic 
field (IWaldl. Il974l : bicak & Dvofakl . ll976b . 

We start out from the F^v describing the test field with asympt otic form of a 
general (i.e. not necessarily parallel) uniform magnetic field given by lBicak & Janis 
(11985) . Due to the axial symmetry of Kerr space-time only two components of 
asymptotic field were considered in that paper without any loss of generality (asymp- 
totic components Bq (parallel) an d B] (equatorial) to be specific). We rewrite com- 
ponents of EM tensor (eq. (A3) of lBicak & Janisl . ll985h denoting Bx = Bi,Bz = Bq 
and splitting the result into two parts according to the asymptotic component. We 
obtain the asymptotically perpendicular part of the field: 



Ff/ =BjcaMrZ~^A~ ^ sin 6 cos 6 [(r^ - 2Mr^ + ra^ ( 1 + sin^ 9) + 2Ma^ cos^ 9) cos \\f 

-a{r^ - 4Mr + ( 1 + sin^ 9) ) sin \|/] , 
F^q' =B^aML~^ (r^ cos 29 + cos^ 9) {a sin \|/ - r cos \|/) , 
F^^^ ^B^aMir ^ sin 9 cos 9 (a cos \|/ + r sin \|/) , 
F^^ =-5_^(<2C0S\|/+rsin\|/) 

-B^aA~^ [(Mr-a^sin^9)cos\|/-a(rsin^9 + Mcos^9)sin\|/] , (2.27) 
F^^ =-Bjc sin 9 cos 9 [(r - Ma^A~ ^ ) cos \|/ - a( 1 + rMA" ^ ) sin \\f] 

+ asm^ QFf^/, 

Fq^ =Bjc [(r^ sin^ 9 + Mr cos 29) cos \|/ - a(r sin^ 9 + Mcos^ 9) sin \\f] 
+ (r^ + a^)B^Mir^ cos 29 + cos^ 9) {a sin \|/ - r cos \|/) 

and the part which approaches uniform field aligned along the axis: 

F^,' =B,aML-^ {?-a^ cos^ 9) ( 1 + cos^ 9) , 

=25,aMrZ-^sin9cos9(r^ -a^), (2.28) 
Fr4 =B,r sin^ 9 + B,a^ sin^ 9MZ"^ {r^-a^ cos^ 9) ( 1 + cos^ 9) , 
Fq^ =5^Asin9cos9 + 2(/ -a'^)5^MrZ"^sin9cos9, 

where we use the azimuthal coordinate \|/ of Kerr ingoing coordinates, which is 
related to Boyer-Lindquist coordinates as follows: 

V = (p + — — In ^— ^ , (2.29) 

r+ — r_ r — r_ 

with r±=M± \JM^ — a^ denoting the outer and the inner horizon. We notice that 
lim^^oo\|/ = (p. 
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As we shall introduce a drift of the black hole in the general direction we lose 
axial symmetry and need to consider all spatial components of the asymptotic mag- 
netic field. We obtain Fj^y (which may only appear due to nonzero drift) by rotating 
Fjjy along the z-axis by angle | - i.e. F^J = F^^y ((p ^ cp — ^ By) which only 
causes sin\|/ — cos\|/ and cos\|/ sin\|/. 

Since the drift shall induce uniform electric field in the asymptotic region we 
need to have appropriate F^v handy. We get them easily by performing dual 



transformation of FjJ'^'\ Dual transformation is carried out as follows: 



1 



2" '^yuvap? 

where £^va|3 is the Levi-Civita tensor whose components are given as: 



(2.30) 



(2.31) 



with [/iVaP] denoting completely antisymmetric symbol. Determinant of the Kerr 
metric is g = 8tt8rr8m8(^(^ - 8lt§rr8QQ = - sin^ 6 Z^. 

Performing the dual transformation we immediately obtain EM tensors with 
desired asymptotics of uniform electric field: 



(2.32) 



In the explicit form we get: 



F^' 



sine E (Fet-(B.,,,, ^ -E,,^,,)^'^^ +Fe^--(B,,,,, ^ -E,,y^,)g'^) 



sin6 Z 



sineZF;e^-(B,,y,,^-£,,,,jrg' 



(2.33) 



sinGSF^/*^ 



sineZ(C-(B,, 
sin9 Z 



-E,^y^^)gV+F^;'-^{B,,y, 



)8")8 



Ex,y,z)8 



Now we are fully equipped to construct any asymptotically uniform test field 
on the Kerr background just by linear superposing of the above EM tensors. As 
we are concerned in constructing F^v which describes the test field around the 
black hole drifting through asymptotically uniform magnetic field in general di- 
rection, we shall employ Lorentz transformation to find the correct asymptotic 
components of such a field. Once obtained we just use them to replace the orig- 
inal "non-drifting" quantit i es Er ,E^,Ey,Br,By and B^. Matrix of general Lorentz 
transformation is (iJacksonl . ll999b : 



( y 

-yvy 



-yv- 

1+^ 



(T-l)vj 



{y-l)VxVy 

{y-l)vxVz 



yvy 

{y-l)VxVy 

(Y-l)vyV; 
„2 



-yvz 

il-l)vxVz 




\ 



1 + 



(2.34) 
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where v = (v^ + Vy + v^) 2 and y = (1 — v^) 2 . 

B B 

Our original field, F^v = +^^v% has a simple asymptotic form in Minkows- 
kian coordinates, 



jriasymptotic i 
^/uv I 



/o 


















5x 



(2.35) 



VO -B, Oj 



To transform the covariant tensor Fjj^, the inverse Lorentz transformation A^, = 
(A^)~^ would be used. But we realize that the Boyer-Lindquist coordinate system 
which we use to perform all the calculations (and also to express the EM tensor of 
the final field) is centered around black hole and the rest frame of the black hole 
is thus our "laboratory" reference frame. As we consider a drift of the black hole 
against the field eq. (12.351) , we need to perform inverse Lorentz transformation. 
Quantities B^c and B^ appearing therein would be primed in standard notation. For 
inverse transformation of covariant tensors we use original A^. Thus for "drifting" 
F^v we have (denoting F^/y' that of eq. (12.351) ): 



77asymptotic 7-. asymptotic a z^' \ v' 



which may be written in the matrix formalism as follows: 



7^ asymptotic 1 



I asymptotic 1 



IIA^II 



I pasymptotic i 



(2.36) 



|A;||, (2.37) 



and results in: 



P asymptotic i 







VyjB^ 



withN=^{v,B, + v,B,). 



VyjBz 


-yB, + v,N 

-VyN 



yB,-v,N 


-yB, + v,N 



1 







(2.38) 



Final step of the derivation is thus substitution of "non-drifting" quantities 

E B 

Ex^y^z and Bxj.z in the tensors F^y'-'^ and F^v by Lorentz transformed values from 
eq. (12.38b and superposing the components to acquire general EM tensor describ- 
ing the field around the Kerr black hole drifting through asymptotically uniform 
magnetic field of general orientation: 



+ /^^ {By ^ -VyN) + F^^ {B, ^ yB, - v,N) . 



-VyyB,) + V (Ey ^ jivxBz - VzB,))+ 
y VyyB,) + F^^ [B^ ^ yB, - v,N)+ 



(2.39) 
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2.4 Choice of the tetrad 



To settle our notation of the tetrad formalism we review basic relations between 
the tetrad basis vectors e^^^^ and dual basis 1 -forms e^^^ 



(a)^ 



r|(^)(v)=^a(3^('^)^fv)' 
.?A^v — 11(a) 5 



(2.40) 

(2.41) 

(2.42) 
(2.43) 



i.e. tetrad indices are raised/lowered using the Minkowskian metric tensor 

Tetrad basis vector e^^^^ is defined by the four- velocity of the observer carrying 
the tetrad ^^^^ = u^. Covariant components are thus related to the tetrad basis as 



follows Uu 



.(0 



2.4.1 Locally non-rotating frame 

One of the standard tetrads which are often used when dealing with Kerr geometry 
is the one carried by the zero angular momentum observer (ZAMO). ZAMO's rest 
frame is usually called locally non-rotating frame (LNRF) since it seemingly 
suppresses intrinsic rotation of the ge ometry. LNR F basis vectors and dual basis 
1 -forms may be expressed as follows dBardeen et al.[ll972b 



4) 



[r) 



^(e) 



,(0 



^ ~ Ai/22:i/2 

Ai/2 
0,^,0,0 

0,0,-^,0 



[i,o,o,a], 



0,0,0, 



Ei/2 



Ai/2sin0 
^1/2^1/2 



Ai/2 

Zl/2 

0,^,0,0 

o,o,2:l/^o 

Ai/2sine 



,0,0,0 



El/2 



[-Q,0,0,1], 



(2.44) 
(2.45) 

(2.46) 

(2.47) 

(2.48) 

(2.49) 
(2.50) 
(2.51) 
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where f2 = — The way in which ZAMO "fits the geometry" may be expressed 

by relation g^^u^zAUO ^ <?^^^zamo which appears useful in the calculations. For 
a detailed discussion of the properties of LNRF and other stationary frames see 
Semerakl([l993 ). 



Elegant way to calculate tetrad basis vectors of any other frame is Lorentz 
transformation of LNRF basis vectors. In order to do this one needs to express 
linear velocity v^^^ of the new frame relative to ZAMO: 

This speed defines the Lorentz boost with factor 7= (1 — v^)~^/^. Change of the 
basis matrix for this transformation is the inverse Lorentz matrix (A^)~^ For the 
transformation of basis 1 -forms the direct A^ is used. 

2.4.2 Frame of the free-falling observer 

We shall calculate the basis vectors of the frame attached to the inertial observer 
who is free-falling from the rest at infinity (FFOFI). Substituting L = and E = 
m (part i cle's rest energy in geometrized units) into Carter's equations of motion 



(ICarted . ll968h vields 



t _^ r _ /2Mr(r2+^y/^ ^ _ ^ _2Mar 

^FFOFI — ^5 WpFOFI — ~ \ ^2 ) ' ^FFOFI — ^5 %FOFI " 



AZ ' 
(2.53) 



which we substitute into eq. (12.521) obtaining single nonzero component v^^^ = 
2r{r +a ) \ r^^^ boost is indeed purely radial: v^^^ = although WppQpj ^ 



^LNRF- Lorentz factor reads Y= J/i^i/i • New (primed) basis is then obtained 
straightforwardly from the LNRF (unprimed) basis by the matrix multiplication 




(2.54) 

We read FFOFI tetrad basis vectors from the above-given matrix equation. Cor- 
responding 1 -forms are obtained from the analogous equation by arranging their 
components in a very same manner as those of basis vectors and switching to di- 
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4) 



,(0 



A /2Mr{r^ + a^ 



E2 



1/2 



IMar 
AZ 



(2Mr(r2 + a2)) 1/2^1/2 ^1/2 2Mflr(2Mr(r2 + a2))i/2 

1 ^1 



' 2: 



ai:aV2 



0,0, 



Zl/2 



,0 



2:1/2 

0,0,0, -yj— - 



i,i^MK!^±^,o,o 



A172 
o,o,2:'/^o 



,0,0 



2:1/2 



a,o,o,i], 



(2.55) 

(2.56) 
(2.57) 
(2.58) 

(2.59) 

(2.60) 
(2.61) 
(2.62) 



we notice that e^g^ and e^'^^^ are common to both LNRF and FFOFI frames. 

Both ZAMO and FFOFI observers are physical everywhere above the horizon. 
In the following, however, we shall set up an observer which is restricted to the 
equatorial plane 9 = 7i/2. 

2.4.3 Frame of Keplerian observer 

As we are primarily interested in astrophysically relevant situations we will employ 
the orthonormal tetrad carried by the inertial observer on the circular Keplerian 
orbit around the black hole (KEP tetrad). Such an orbit is specified by the values 
of constants of motion - by specific angul ar momentum L = u,n and specific energy 
E = —u, which are expressed as follows ( Bardeen et aI.LIl972 ): 



-2Mr±a\/Mr 



r\/ — 3Mr ± la\fMr 



r{r'^ -3Mr±2a\/Mr) 



(2.63) 



where the upper signs are valid for the prograde (direct) orbits and the lower ones 
for the retrograde (counter-revolving) orbits. Such a tetrad is physical only above 
marginally stable orbit r^as which represents a radial boundary for the stationary 
geodesic motion in the equatorial plane: 



:M(^3 + Z2T V(3-Zi)(3 + Zi + 2Z2)j, 



(2.64) 
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where Z 



1 



A/2 



1/3 



1/3' 



and Z2 = 



3^!i -L 72 pn- 
^ -t-z,j. ro 



sition of the marginally stable orbit is plotted as a function of spin a in fig. IZT 




Figure 2. 1 : Position of the marginally stable orbit for both co-rotating and counter- 
rotating equatorial circular geodesies as well as the position of the outer horizon is 
plotted as a function of BH spin parameter a. 

Angular velocity of a circular orbit is i^Kep = M-^/V/^±a ^^^^^ upper 
signs are for prograde orbits and the lower ones for the retrograde orbits. Lin- 
ear velocity eq. (12.52b of the orbiting tetrad as measured by ZAMO observer in 
LNRF reads v^f) = (^Kep -^^)- Lorentz factor of this azimuthal boost is 
Y = (1 — [v^'P)]^)"'/^. By means of Lor entz transformation of L NRF basis we ob- 
tain basis vectors of the orbiting frame dYokosawa & Inui . 2005 '): 



/ 4 \ 1/2 



(r) 



1/2 



[0,1,0,0], 



(2.65) 



1/2 



,0,0, 



sin 



+ v 



AL 



1/2- 



where we define A = (r^ + a^)^ — a^Asin^ to express the coordinate angular ve- 
locity of LNRF (i.e. angular velocity of the frame-dragging) Q = ^Mr. 

For r < Tins there are no more circular orbits. Thus we suppose that the orbiting 
observer who reaches this limit performs a free fall to the black hole keeping the 
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values of the constants of motion corresponding to the marginally stable orbit at 
(FFO tetrad). He is falling with = ^(^ms) and L^s = ^(^ms) given by eqs. 
(12.631) - (12.641) . Having fixed Ut{r < r^s) -Ems, W(p(r < r^s) ^ms and we = 
we get radial component Ur easily from the normalisation condition u^Ujj = — 1. 
Contravariant components of the 4- velocity are then: 

= ^ M^^ + + ^Ma^)Ems - 2MaLms] , 
m 



1 



u'- = -^^ [r(r2 + a2) + 2Ma^]El^ - AMaEr^,L^, - (r - 2M)Ll^^ - rA, 

(2.66) 

= 0, 

m'P = ^ [2M2fl£^s + M{r - 2M)L^,] . 

Spatial 1 -forms of the tetrad of this falling observer may be expressed as follows 
(lDovciakl . E004lv 

(r) ^ 

^==^=K[w^,w^,0,W(p] + [0,1,0,0]), 



[0,0,r,0], (2.67) 
^"^,0,0,^^. 



1 + U^Ur 



We conclude the we are now equipped with three distinct tetrads to be applied 
in the equatorial plain, namely non-inertial LNRF and two inertial tetrads FFOFI 
and KEP+FFO. Outside the equatorial plain only LNRF and FFOFI may be ap- 
plied. 



2.5 Structure of the electromagnetic field 
2.5.1 Stationary electromagnetic field 

First we shall review several aspects of the stationary, i.e. non-drifting, asymptoti- 
cally homogeneous test fields given by eqs. (12.271) and (12.281) . Since the stationary 
fields of this type has already been discussed thoroughly in the literature (see e.g. 



Bicak et al.Lll989l , and references therein) we will focus mainly on the comparison 
of the field line structures resulting from the alternative definitions of the field com- 
ponents which were given in sec. 12.21 We also discuss the choice of the test charge 
four-velocity profile (ZAMO versus FFOFI basically). Besides that we shall in- 
troduce various techniques to visualize the vector fields in both, two-dimensional 
plane sections as well as the stereometric projections of three-dimensional space. 



CHAPTER 2. MOTION IN BLACK HOLE MAGNETOSPHERES 



24 



Expulsion of the aligned magnetic field (Meissner effect) 



We briefly revisit the issue of the expulsion of the axisymmetric stationary mag- 
netic field out of the horizon of the extreme Kerr black hole (a = M) which is known 
as Meissner or Meissner-type effect in the analogy with similar effe ct which super- 
conducting bodies exhib i t upon the external magnetic fields (e.g. IPovciak et al.L 
200d : iBicak & Ledvinkal . l200d and references therein). In particular we shall dis- 
cuss the role of the definition of the lines of force and the observer dependence of 
this effect. 




Figure 2.2: Meissner effect captured by AMO components which translate the van- 
ishing of the F0(p at the horizon directly into the radial component of the magnetic 
field ^AMO- Upper left panel shows the case a = (Schwarzschild limit) where the 
field lines penetrate the horizon completely indifferently. Increasing the spin the 
lines of force start to bend slightly which becomes obvious for higher spin values 
only. Upper right panel captures the situation for a = 0.9 M. For a = 0.99 M the 
effect becomes more apparent (bottom left panel) and in the extreme case a^M 
the magnetic field is expelled from the horizon completely (bottom right panel). 



Scalar magnetic/electric flux ^^j^ through the given surface S may be calcu- 
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lated as another quantity revealing the structure of the field. In particular the gen- 
eral expression for the magnetic flux 



(2.68) 



reduces to the simple form (^j^ = / ^69 d9dcp if we choose the surface to be a part 
of the sphere of constant r. If we consider the case of aligned field of the strength 
given by eq. (12.281) and further specify the surface to be an axisymmetric polar 
cap determined by the polar angle 6 we may integrate straightforwardly to obtain 
(l)ni = TiBz [A+ ^^(r^ — a^)] sin^9. The flux through the horizon vanishes in the 
case of extreme black hole and magnetic field is fully expelled fro m the horizon. 



Poloidal section of the surfaces of the constant flux discussed by iDovciak et al. 



(l200d ) coincide with the AMO field lines since they both directly reflect the be- 
haviour of the Fq^ which is zero at the horizon of the extreme BH. Expulsion of 
the magnetic field in AMO definition is visualized in figs. 12.21 and [231 The latter 
employs the linear integral convolution (LIC) method (Shambo , 2005) which en- 
codes the field structure into the texture resembling iron filings. The strength of 
the field is expressed by the color scale. Visualisation of the vector fields using 
LIC method proves especially useful when dealing with complex field structures 
for which the usual field lines do not keep up with the sudden spatial changes (e.g. 
dense zigzag structures, crack propagation etc). 




-1 I 



Figure 2.3: Alternative vis ualization of th e field lines employing linear integral 
convolution (LIC) method (IShambd . l2005h . Strength of the field is encoded by a 
colormap: brighter color means stronger field. Amplification of the magnetic field 
in the equatorial zone in the vicinity of the horizon is observed. 

In figs. 12.41 and [231 we review behaviour of the magnetic field lines of various 
definitions for both ZAMO and FFOFI velocity profiles of the test charge in the 
case of extreme BH a = M. We observe that coordinate components B^ for both 
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ZAMO and FFOFI four- velocities exhibit the Meissner effect and magnetic field is 
expelled from the horizon (top panels of figs. I2.4l and [231) . 

Large-scale bending of the field lines is necessarily artificial since both ZAMO 
and FFOFI are asymptotically static and the distant test charge thus shall perceive 
asymptotic shape of the field which is aligned uniform magnetic field in this case. 
Bending is caused by the lack of normalization of the coordinate basis vectors 
which led us to the definition of the physical components given by eqs. (12.141) - 
(12.161) . In the second rows of figs. 12.41 and [23] we observe, however, that physical 
components do not exhibit the expulsion of the magnetic field; neither for ZAMO 
nor for FFOFI. This is caused by the fact that the definition of physical compo- 
nents incorporates the coordinate singularity at the horizon directly into the field 
component ^physical commented above. By exploring the behaviour of the un- 
derlying expressions we learn that in ZAMO case the radial component has finite 
limit at the horizon lim^^^^ ^physical ^ ~ (f+cos^e)^ while in FFOFI case it diverges 
as ^physical ^ We thus coufirm that physical components ^l^y,^^^^, physical 
are inappropriate when working in Boyer-Lindquist coordinates. 

Renormalized components which we suggested to use to overcome this prob- 
lem actually eliminate the divergence of the radial component at the horizon (in the 
limit r ^ r+ renormalized components behave in the same manner as the coordi- 
nate ones) and also provide the proper asymptotic behaviour (see the third rows in 
figs. I2.4l and [23]) . Finally we notice that FFOFI observer in his frame observes the 
expulsion (bottom panels of fig. 12.51) while the ZAMO observer sees the magnetic 
field lines penetrating the horizon (fig. 12.41) . Expression for 5^^^ involves factor 
terms el ^u^ — ^u^ and ^u^ — .u^ which are both zero in the limit r ^ r+ in 

(r) (r) (r) (r) + 

the FFOFI case and B^^^ thus cancels out. For ZAMO which has — 0, how- 
ever, these terms do not cancel and thus allow B^^^ to acquire nonzero value at the 
horizon of the extreme Kerr BH. 

On the other hand 7^ causes that FFOFI measures nonzero azimuthal com- 
ponent of the magnetic field given as ^ppopi = —eJ^^^F^^rU^ which has real valued 

limit lim^^^^ ^ffofi ^ ~ (1 +00^29^/2 horizon of the extreme Kerr BH. We fur- 
ther notice that renormalized components of the field felt by the FFOFI test charge 
encompass the azimuthal component ^renormalized ^ -w'^rsinG (g^^ *F(pr + */^rr) 
which actually diverges as ^renormalized ^ [r-r+y r ^ r+ in the extreme Kerr 
geometry a^M where r+ = M. 

Due to the presence of the azimuthal component the FFOFI field lines are not 
confined to the given poloidal plane (p = {const ^ const + Ti} and the ordinary two- 
dimensional plane plots of the lines of force thus represent a projection of the 
vector field rather then the true field lines. Under these circumstances one should 
occasionally check the stereometric projection of the field lines in order to reveal its 
azimuthal component since it may change the picture profoundly. We emphasize, 
however, that the axial symmetry is maintained and the azimuthal component is 
independent of (p. 
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Figure 2.4: Magnetic field felt by ZAMO test charge orbiting above the horizon 
of the extremal Kerr black hole immersed into the aligned field. We compare (top 
to bottom) coordinate, physical, renormalized and tetrad components of the field 
showing both, the overall shape of the field lines as well as the detailed bahaviour 
in the vicinity of the horizon. 
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Figure 2.5: Magnetic field felt by the FFOFI test charge freely falling onto the 
extremal Kerr black hole immersed into the aligned field. We compare (top to bot- 
tom) coordinate, physical, renormalized and tetrad components of the field show- 
ing both, the overall shape of the field lines as well as the detailed bahaviour in the 
vicinity of the horizon. 
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We introduce stereometric plots by showing the simple AMO case (which has 
no azimuthal component since F^q = 0) for the two values of spin a = 0.9M and 
a = M in the upper panels of fig. 12.61 These may be directly compared with the 
corresponding plane figures of fig. 12.21 Then we move to the FFOFI case where the 
azimuthal component is present. In bottom panels of fig. 12.61 we compare FFOFI 
tetrad components with renormalized FFOFI components. The latter exhibit diver- 
gence in the azimuthal component at the horizon which causes strong winding of 
nearby field lines while the tetrad field lines are twisted only partially, shifted by 
the finite angle in the r ^ r+ asymptotics. 

In the following sections we shall preferently work with the tetrad components 
of the fields as they provide a consistent interpretation of quantities measurable 
by a physical observer. When occasionally switching to the coordinate basis we 
use renormalized components for their properties. Sometimes it appears useful to 
visualize the components of the EM tensor itselves without notion of the specific 
four- velocity of the test charge - then the AMO components are employed. 
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Figure 2.6: Meissner-type expulsion of the magnetic field from the horizon of the 
extreme Kerr black hole. In the upper left panel for a = 0.9M we observe slightly 
bent AMO field lines piercing the horizon. Upper right panel shows the extreme 
case a = M for which the AMO magnetic field is fully expelled. Bottom panels 
show the field as felt by FFOFI. In the bottom left panel we observe tetrad compo- 
nents of the magnetic field with nonzero azimuthal component which approaches 
finite nonzero value on the horizon. Renormalized FFOFI components, however, 
exhibit divergence in the azimuthal component as we approach the horizon. In 
the close neighborhood of the horizon we observe strong azimuthal winding of the 
field lines (bottom right panel). Width of the "tubes of force" is chosen arbitrarily 
to optimize legibility of the plots. 
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Non-aligned uniform magnetic field 

In the following we shall in brief consider also a field component perpendicular 
to the rotation axis. We therefore lose the axial symmetry. Lines of force are not 
anymore confined to the poloidal sections as they were in the case of stationary 
observers/test charges in the axisymmetric situation. We therefore plot the sections 
in couples of (x,z) and (y^z) plane projections (by setting (p = {0,7i} and (p = 
{71/2,371/2}, respectively) of the field lines to capture at least fractionally their 
azimuthal dependence. Stereometric projections are often employed to provide 
more insightful view on the complex field structure. 




Figure 2.7: AMO components of the asymptotically homogeneous magnetic field 
perpendicular to the symmetry axis (B^ = 0). Left panel depicts the Schwarzschild 
limit <3 = for which the black hole does not affect the field at all as we have 
already observed in the analogical aligned case in fig. 12.21 The behaviour of the 
field lines, however, differs profoundly once the rotation is switched on. In the 
right panel we observe the field lines twisted by the extreme spin a = M. Unlike 
the aligned case, here we observe no expulsion of the field. 

First we mention that by eq. (12.271) the component does not diminish at 

the horizon of the extreme Kerr black hole as Fq^ does. Magnetic field lines thus 
penetrate the horizon; effect of the magnetic expulsion is not present in the AMO 
field components which we illustrate in fig. 12.71 where we compare AMO mag- 
netic field lines for two extremal values of spin, namely a^O and a = M, setting 
purely perpendicular magnetic field on the background: B^ = B^ = 0. In the 
Schwarzschild limit the field lines are not affected by the black hole since the F^y" 
components reduce to its asymptotic form immediately by setting a = 0. Setting 
nonzero spin the symmetry is lost and the field lines are twisted in a rather compli- 
cated way as we observe in the right panel of the fig. 12.71 In general the field lines 
may penetrate the horizon even in the extreme case a — M which we confirm also 
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for other definitions of the field lines and for both ZAMO and FFOFI four- velocity 
profiles by exploring the underlying formulae. 

We further notice that in the equatorial plane the asymptotically perpendicular 
magnetic field (Bjc ^ 0, 5^ = 0) will have zero longitudinal component = = 
regardless the field definition also in the vicinity of the black hole (only equatori- 
ally vanishing terms from (12.271) appear in the definition of B^ in eq. 12.111) . In other 
words the field lines reside in the equatorial plane in this special case. This might 
not be that visually apparent from the right panel of fig. 12.71 Therefore we add a 
series of equatorial sections in fig. 12.81 to reveal the field in the equatorial plane. 
Three equatorial observers, namely FFOFI, co-rotating KEP+FFO and counter- 
rotating KEP+FFO are involved in the comparison. We observe that free fall with 
L = (FFOFI) results in the aligned magnetic field lines in the Schwarzschild limit 
which we observed also in AMO components in fig. 12.71 as a result of F^y taking 
its asymptotic form for a = 0. On the other hand, nonzero angular momentum 
of KEP+FFO observers which are freely falling with L = L(rins) under r^^ causes 
twisting of the field lines even in the Schwarzschild case since the four-velocity 
and other tetrad vectors components do not take their asymptotic form with a = 
although the F^v does. When the spin is employed the lines of force twist in a rather 
complex way regardless the choice of the observer. Complex layering of the field 
lines emerges in the narrow zone above the horizon. In order to stretch this tiny 
yet very interesting area and make it more convenient for visualisation we rescale 
the radial component introducing the dimension-less coordinate R = which 
squeezes the horizon into a single point R = and radial infinity into the circle 
7^ = 1. Rescaled coordinate R is employed in bottom panels of fig. 12.81 

In the series of stereometric projections in fig. l2.9l we observe the interplay be- 
tween the aligned and perpendicular magnetic field components using AMO def- 
inition of the field on the extreme Kerr background. We start with the aligned 
situation where the expulsion of the field is observed and gradually increase the 
perpendicular Bjc component to see that it generally causes intricate deformations 
of the field lines and allows them to penetrate the horizon even in the extreme Ken- 
case. 

In fig. 12.81 we observe that the field twisting caused by the interplay between 
the black hole rotation and the horizontal component of the prescribed magnetic 
field is accompanied by intensive layering of the magnetic field in the narrow zone 
above the horizon. In fig. 12. 101 we compare the detailed structure of the magnetic 
layer zone as measured by several distinct observers. We observe that magnetic 
layers are present regardless the choice of the frame although their actual shape 
changes as we switch between them. It is thus apparent that magnetic layers do not 
emerge as a observer effect. They rather originate directly from the interplay be- 
tween the frame-dragging effect and prescribed magnetic field which is confirmed 
by observing layered magnetic structures also in observer independent AMO com- 
ponents. 

Simulations of non-vacuum magnetospheres carried out in the framework of 
force-free electrodynamics (FEE) often reveal similar magnetic layers in which the 
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Figure 2.8: Equatorial behaviour of the magnetic field with purely perpendicular 
asymptotics (Bjc 7^ 0, 5^ = 0) for three slightly distinct tetrads. We remind that 
FFOFI stands for the observer which is freely falling from the rest at infinity (an- 
gular momentum L = 0) while KEP+FFO is the tetrad attached to the Keplerian 
observer being stable above the marginally stable orbit r^s (eq. 12.641) and falling 
freely below this orbit keeping L{r^s). The top row presents the FFOFI tetrad, 
the second the prograde KEP+FFO and the third is for the retrograde KEP+FFO. 
Schwarzschild limit a = is shown in the first column, middle column corresponds 
with a = 0.9M and the last one represents the extreme case a^M. In the bottom 
row the rescaled dimensionless radial coordinate R = is used to stretch the re- 
gion close to the horizon which is of the utmost interest. Panels of the bottom row 
show the fields measured by FFOFI and the both, co-rotating and counter-rotating 
KEP+FFO for extremal BH case. 
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Figure 2.9: Series of stereometric projections illustrates the impact of the perpen- 
dicular component B^c upon the originally axisymmetric magnetic AMO field lines 
of Bz origin in the case of extreme Kerr BH. Upper left panel shows aligned case 
(Bjc = 0) where the field is expelled, upper right ^ = 0.1, bottom left ^ = 0.5 and 
bottom right with Bx = B^. Perpendicular component Bx generally causes complex 
twisting of the field lines and allows the field to penetrate the horizon even in the 
extreme case. 
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Figure 2. 10: Magnetic layers above the horizon of extreme Kerr BH that develop as 
a consequence of an interplay between the spin induced frame-dragging effect and 
perpendicular component of the prescribed magnetic field. Fraction of the layered 
zone is shown. We remind that in this case the field lines reside in the equatorial 
plane. Panels show the field measured by (from the top left to the bottom right): 
FFOFI observer, co-rotating Keplerian observer, observer freely falling from the 
counter-rotating marginally stable orbit and ZAMO. Last two panels show observer 
independent AMO components of the field in which the layering is also apparent 
proving it to be an intrinsic feature of the system rather than a pure observer effect. 

direction of the field changes sharply. FFE assumes that plasma is streaming solely 
along t he magnetic field lines and such layers thus become current sheets. See 
Spitkovsky, (.2006 ) for the simulation of pulsar magnetosphere in which the current 
sheets develop. 

Electric field 

Until now we have paid hardly any attention to the electric fields felt by the test 
charges in a given setup of Kerr black hole embedded in the asymptotically uni- 
form magnetic field with general orientation with respect to the rotation axis. To 
correct this we present a summary of nonzero electric field components for var- 
ious field definitions, observers and values of the spin parameter a in tab. 12.5.11 

B B 

We distinguish whether F^y or F^J components (or actually both) act as a source 
terms of a given electric field component. We treat separately the case of equatorial 
plane 6 = f where part of the components diminish. Exploring tab. 12.5. II we notice 
particularly that nonzero electric field is measured even in the Schwarzschild limit 
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Figure 2.11: AMO components of the electric field above the extreme Kerr BH. In 
the left panel we depict the aligned case = B^c^O for which no azimuthal 
field components is measured. In the equatorial plane this field becomes purely 
radial. On the other hand the non-aligned situation B^^Bjc^ presented in the 
right panel brings azimuthal electric component into the play. 

Stereometric projections of AMO electric fields above extreme Kerr BH are 
presented in fig. 12.111 to illustrate the profound changes in the field line structure 
which accompany the onset of perpendicular magnetic field component B^ which 
disrupts the symmetry of aligned case. 

In the series of plane sections in figs. 12.121 and [2.13l we compare the structures 
of the electric fields as measured by ZAMO and FFOFI observers. First we present 
the situation of aligned field (B^ = 0) in which all poloidal sections of the fields 
coincide due to the symmetry (top panels of fig. l2.12l) . We notice that in the FFOFI 
tetrad (unlike the ZAMO tetrad) the electric field is expelled from the horizon of 
the extreme BH. Once the non-aligned field is considered (by setting Bjc ^ 0) the 
symmetry is lost and poloidal sections for different values of cp generally differ 
(see the bottom panels of fig. l2.12l and fig. l2.13l) . The effect of electric expulsion is 
gone once the non-aligned field is considered. 

Stereometric projections in fig. l2.14l reveal that in the renormalized components 
the expulsion of the electric field out of the horizon of the extreme Kerr is present 
not only for the FFOFI test charge but also for the test charge with ZAMO four- 
velocity (while in the ZAMO tetrad the field lines are not expelled). Azimuthal 
component which is present in the FFOFI case causes winding of the field lines 
while the ZAMO field lines remain confined to a given poloidal sections. Bottom 
panels of fig. l2.14l illustrate that setting B^^O generally disrupts the symmetry and 
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Table 2.1: Summary of the electric field components depending on the definition 
of the field, spin of the black hole and employed preset magnetic field (which may 
be aligned or not). We treat separately the situation in the equatorial plane 6 = | 
since the field tends to simplify here profoundly. Separate columns labeled and 
Bjc suggest for which components of measured electric field is given component 
of the asymptotic magnetic field being "responsible". In other words electric field 
components found in e.g. B^ column are "generated" by Bx. We notice particularly 
that in the case of free-falling test particle (for both coordinate and tetrad compo- 
nents) as well as for Keplerian four- velocity we measure nonzero components of 
electric field even if there is zero spin (Schwarzschild limit). We emphasize that the 
fact that we arrive at the very same components being nonzero (though not of the 
same values!) both in the coordinate frame and the tetrad frame for a given four- 
velocity of the test charge is not in any sense automatic or generally guaranteed as 
we could easily think of such F^v for which it would not be the case. 
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allows the field lines to penetrate the horizon regardless the spin value. 

The expulsion of the electric field as observed in the FFOFI tetrad is further 
studied in fig. 12.151 where we provide series of poloidal plane sections for the 
aligned background field which differ in the spin value. Such discussion allows 
us to visually reveal the mechanism of the expulsion which takes place in the case 
of the extreme spin value. LIC patterns are also employed to help to reveal this 
effect. 

Electric field is generally very sensitive to the value of spin parameter which 
is not that surprising since the field itself has gravitomagnetic origin by which we 
mean that it arises from the interplay between the background magnetic field and 
spin induced frame-dragging effect. As a consequence we measure no AMO elec- 
tric field in the Schwarzschild limit since all Fti vanish with a = 0. Naturally this 
does not necessarily mean that all observers will measure zero electric field. We 
have seen, for instance, that even in the case of zero spin the non- stationary FFOFI 
observer may observe nonzero electric field induced due to his radial velocity com- 
ponent u^. Nevertheless these electric components asymptotically diminish since 
the observer is asymptotically static. 

The question of the expulsion of the electric field has been addressed in the case 
of aligned background field. Unlike the magnetic expulsion where the radial AMO 

component ^amo ^ r^sinQ ^^^^i^hes at the horizon for a = M here the analogical 
component ^'amo ^ ^^^^ vanish. Thus the electric expulsion is merely an 
observer effect. We have seen that in the ZAMO tetrad the electric field was not 
expelled while in the FFOFI one it was. We recall that also the magnetic field is 
fully expelled in the FFOFI frame. 
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Figure 2.12: Projections of the electric field onto the (x,z), (y^z) and (x^y) planes 
in the case of the extreme Kerr BH with the asymptotic magnetic field. Top six 
panels capture the situation in the aligned setup, = and = M~^. Due to the 
axial symmetry the poloidal projections (x,z) and {y^z) coincide in this case. We 
compare ZAMO tetrad components in the first row with the FFOFI tetrad compo- 
nents in the second one. FFOFI differs from ZAMO most strikingly in measuring 
nonzero azimuthal components E^^\ Second series of six panels represent the sit- 
uation of transversal orientation of the magnetic field, = and B^ = 0. Axial 
symmetry is lost, projections onto different poloidal planes generally differ. We 
compare ZAMO tetrad components in the third row with the FFOFI tetrad compo- 
nents in the fourth one. 
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Figure 2.13: Projections of the electric field onto the (x,z), {y,z) and {x,y) planes 
in the case of the extreme Kerr BH with the asymptotic magnetic field components 
= and = M~^. We compare ZAMO tetrad components in the upper 
panels with the FFOFI tetrad components in the bottom ones. ZAMO measures 
azimuthal field E^^'^ outside the equatorial plane only, see tab. 12.5. II for the sum- 
mary. 
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Figure 2.14: Electric field is expelled from the horizon of the extreme Kerr black 
hole in both ZAMO and FFOFI renormalized components in the case of aligned 
magnetic field on the background (B^ = 0). In the case of ZAMO (upper left 
panel) there is no azimuthal component while the FFOFI observer measures also 
azimuthal electric field causing winding of the field lines (upper right panel). How- 
ever, considering B^y^O brings azimuthal electric component also for ZAMO as 
we observe in the bottom panels. The bottom left shows the situation for |^ = 0.1 

while the right one for ^ = 0.3. Electric field is not expelled in general, some field 
lines penetrate the horizon. 
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Figure 2.15: Evolution of the FFOFI measured electric field in the case of aligned 
background magnetic field (B^c = 0) when the spin is gradually increased is shown 
in the series of poloidal plane sections. We begin with the low value a = O.IM 
(upper left panel) which is increased to the medium value a = 0.5M (upper right 
panel), high value a = 0.9M (middle left) and eventually to the extreme value 
a = M in the middle right panel. We observe the effect of the expulsion of the 
electric field in the extreme case. Bottom panels show LIC patterns reflecting the 
expulsion. Bottom left panel shows a = 0.9M case with only partial expulsion, the 
right one captures the complete expulsion occurring for a^M. Color bar suggests 
the strength scale of the field (units are arbitrary). 
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2.5.2 Electromagnetic field around drifting black hole 

In this section we shall discuss the electromagnetic fields around Kerr black hole 
which is drifting in a general direction through the oblique uniform magnetic field. 
Components of EM tensor F^v describing such field were derived in sec. 12.31 in 
terms of Lorentz transformation of former non-drifting solution given by eqs. I2.27[ 

Resulting system provides highly simplified though not unrealistic model of 
conditions occurring in a diluted gaseous medium in the vicinity of the inner edge 
of the accretion disk. This region may be threaded by the external large-scale 
magnetic field which we model by oblique uniform field in our setup. Gravitomag- 
netic effect acts on the field enriching profoundly its structure as we have seen in 
previous section. Observer dependence of the EM field components has also been 
discussed previously. In addition to these effects we introduce the translational uni- 
form motion of the black hole itself bringing extra parameters to the system. See 
Karas & Kopacek (l2009h for a discussion of astrophysical relevance of our model. 



Lorentz transformation of the field changes its asymptotical form - in gen- 
eral it changes the direction of the former magnetic field and induces electric field 
which was not present in the asymptotic region of non-drifting BH (for asymp- 
totically static observers at least). In the previous section we observed that inter- 
action between BH's rotation and perpendicular component of the magnetic field 
causes complex twisting of the field lines and that a narrow zone of magnetic layers 
emerges above the horizon regardless the choice of the observer. 

In fig. 12.161 we introduce the drift in the case of extreme BH embedded in the 
aligned field in order to detect the drift impact upon the field most clearly. We 
observe that the Meissner effect is suppressed by the drift. We note that with the 
drift velocity increasing the shock front develops in which the field has complex 
layered structure which gradually enhances as the velocity rises. Besides the mag- 
netic layers we newly observe the pair of magnetic null points which emerge for 
the sufficiently rapid drift. Sites of zero magnetic intensity leave the charged parti- 
cles prone to the acceleration by the electric field which is generally not vanishing 
here. The layered field structures whic h surround the neutral p oint present the 
sites of possible magnetic reconnection (iKaras & Kopacekl . 1200 9'. and references 



therein). Magnetic layers and neutral points are further explored in fig. l2.17l using 
Lie patterns with the strength of the field being encoded by the color scale. 

From fig. 12.101 we recall that magnetic layers emerge also in the non-drifting 
case as a consequence of a perpendicular Bjc component of the field. Layering 
accompanied with the drift through the aligned field thus may be attributed to the 
tilt arising from the Lorentz transformation. The neutral points, however, only 
appear in the drifting case. 

In stereometric projections fig. 12.181 we compare magnetic field measured in 
the FFOFI frame to that captured by the ZAMO renormalized components. In both 
the cases the aligned non-drifting field above the extremal BH is fully expelled, 
the only difference being that the FFOFI field is slightly azimuthally twisted (see 
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fig. 12.61) . When the drift is introduced, however, the behaviour of the field lines 
differs strikingly as they are dragged in the opposite direction. FFOFI lines stretch 
in the drift direction in an intuitive way while in the ZAMO renormalized compo- 
nents they are dragged backwards manifesting profound differences between the 
two field definitions. 

Electric field also reacts markedly when the drift is introduced. In fig. 12.191 
we explore its impact upon the FFOFI electric components above the horizon of 
the extremal BH in the poloidal (x,z) plane section. Gradually increasing the drift 
velocity we observe how the topology of the field evolves. We note that without 
the drift the electric field asymptotically decays whereas the drift induces non- 
vanishing electric field which overruns the original one. Thus the global impact of 
the drift is more dramatic compared to what we observed in the case of magnetic 
field. Besides overall changes of the field we also notice progressive layering in the 
narrow region above the horizon. Similarly to their magnetic analogue, the electric 
layers are also present regardless of the choice of the observer and their origin may 
be attributed to the tilting of the originally aligned magnetic field due to the drift. 

We also explore the electric field in the equatorial plane in fig. 12.201 Since we 
start out from the aligned magnetic field and the drift is restricted to be perpen- 
dicular to the axis (v^ = 0) the field lines will reside in the plane (see tab. 12.5.11) . 
The rescaled radial coordinate R = ^-^^ is employed. We compare the tetrad com- 
ponents of ZAMO, FFOFI and both co-rotating and counter-rotating KEP+FFO. 
The drift changes the structure of the field completely including the asymptotic re- 
gion. Although the fields differ from each other noticeably for different frames, in 
all the cases the neutral points are formed (at different location, though) once the 
sufficiently rapid drift is introduced. 

We conclude that the translational motion of the black hole through the uni- 
form magnetic field enriches profoundly the structure of the resulting electric and 
magnetic fields. We observed that in the case of axisymmetric background the drift 
causes a formation of a narrow zone in which the fields (both electric and mag- 
netic) are heavily layered and the layering progressively enhances with the rising 
drift velocity. However, we observe such effect also for nondrifting oblique fields 
(see fig. 12.101) and therefore we attribute the formation of the layers to the perpen- 
dicular component which arises from the field's tilt in the case of drift through the 
aligned magnetic field. Nevertheless only in the drifting scenario we can observe 
such topology of layers which leads to the formation of neutral points. These were 
observed for both, magnetic as well as electric fields. 
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Figure 2.16: In the series of first five panels we observe the impact which the 
gradually increased velocity V;c = 0,0.2,0.5,0.8 and 0.99 of the drift along the 
positive direction of the horizontal axis has upon the FFOFI measured magnetic 
field. Spin is set to the extreme value a = M. Prescribed aligned magnetic field 
is originally expelled from the horizon. When the drift is introduced, however, 
the field lines are allowed to penetrate the horizon. Increasing the drift speed we 
observe that the layered bumper zone develops being further enhanced with rising 
drift velocity. Right panel of the middle row and the left one on the bottom row 
show the magnetic layer structures which arise for = 0.5 in detail. Layered 
patterns further enhance when the speed increases - next panel (middle bottom) 
shows the rapidly moving case = 0.99. Last panel reveals the null point of the 
magnetic field which develops with high drift velocity. 



CHAPTER 2. MOTION IN BLACK HOLE MAGNETOSPHERES 



46 




0.^ 4- n.n O.G 0.7 OA 1 -0.67 -0.64 -0.62 -0.59 -0.57 -0.55 -0.5 



X [M] X [Ml 

Figure 2.17: Drift induced effects upon the structure of the FFOFI tetrad com- 
ponents of the aligned magnetic field in the case of extreme spin value a = M. 
Translational motion is restricted to be parallel with the horizontal axis, ^ 0, 
= = 0. Upper figures reveal the situation for V;^ = 0.5. We observe that a 
narrow buffer zone develops above the horizon in which the field lines have com- 
plex layered structure. Zones of antiparallely oriented magnetic field are brought 
to the close contact here. Bottom panels show the case of extremely rapid mo- 
tion = 0.99. Null point of the magnetic field appears in this case (bottom left 
panel). In the right panel we observe self-similar tightly folded layered structures 
which are considerably enhanced when compared to the case of slower motion in 
the upper right panel. 
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Figure 2.18: Drift induced deformations of the originally aligned magnetic field 
lines above the horizon of the extreme Kerr black hole. FFOFI tetrad components 
in the upper row are compared to the renormalized ZAMO components in the bot- 
tom row. In both the cases the field is expelled for the zero drift velocity (see 
fig. 12.61) the only difference being that FFOFI field is slightly twisted in the az- 
imuthal direction. Drift direction is set to coincide with the horizontal axis of the 
stereometric projection, Vjc > 0, v^; = —Vjc and = 0. Velocity is increased in the 
following sequence: V;c = 0.1, 0.3, 0.7. FFOFI tetrad and ZAMO renormalized 
components exhibit completely different behaviour when the drift is introduced as 
ZAMO field is dragged in the opposite direction compared to FFOFI. 
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Figure 2.19: Poloidal (x,z) plane sections of FFOFI measured electric field 
around drifting extreme Kerr BH. First six panels show the overall changes of 
the field topology as the drift velocity increases in a following sequence v^^ = 
0, 0.1, 0.2, 0.3, 0.5 and 0.99. Bottom series of six panels presents detailed view on 
the progressive layering which takes place in a narrow region above the horizon. 
We observe that as the drift velocity rises the layering enhances profoundly. 
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Figure 2.20: Equatorial behaviour of the electric field around extremal Kerr BH 
drifting in the aligned magnetic field. Employing the rescaled radial coordinate 
R = ^-yi makes the horizon shrink into the single point residing at the origin. First 
column presents non-drifting case, in the second we set v^^ = 0.5, = and in 
the third = 0.5, v^; = 0.5. Four distinct frames are compared in the rows (top to 
bottom): ZAMO, FFOFI, co-rotating KEP+FFO and counter-rotating KEP+FFO. 
We observe that in all the cases the neutral electric points develop as the drift is 
introduced. We stress that as the original field is aligned (B^ = 0) we measure 
= in the equatorial plane provided that = 0. In other words above plots 
present true field lines rather than mere section. 
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2.6 Motion of charged test particles 



Off-equatorial, energetically bound motion of charged particles in strong gravi- 
tational and electromagnetic fields is pertinent to the descriptio n of accretion disk 
coronae around black holes and compact stars. In recent papers (IKovaf et alJ . l2008l . 
2010h we discussed the existence of energetically-bound stable orbits of charged 
particles occurring outside th e equatorial plane, extending thus a large variety of 



iPrasannal. 



complementary studies (e.g.. iBicak et all 119891: IStuchlik et al 
1980; IPrasanna & Senguptaf . 1 1994 : lAUev & Ozdemii , 20021 . and further references 



cited therein). Particles on off-equatorial stable trajectories form a coronal flow that 
is possible at certain radii and for certain combinations of the model parameters, 
namely, the specific charge of the particles, the conserved energy and the angu- 
lar momentum of the particle motion, the strength and orientation of the magnetic 
field, and the spin of the central body. 

We assume that t he magnetic field permeating t he co rona has a large-scale 
(ordered) component (IBisnovatyi-Kogan & Lovelacel . l2007b . In this case, charged 
particles can be trapped in toroidal regions, extending symmetrically above and 
below the equatorial plane and forming two halo lobes. H owever, this trapping 
happens only for certain combinations of model parameters (IKovaf et all 12010) . 

We consider two types of the model setup: a rotating (Kerr) black hole in 
an asymptotically uniform magnetic field paralle l to the symmetry axis (IWaldl . 
19741 : iTomimatsu & Takahashil . l200ll : iKoidj . l2004l : iKoide et al.[ l2006 h. and a non- 



rotating star (desc ribed by the Schwarzschild metric ) endowed with a rotating mag- 
netic dipole field (IPettersonl . 1 1 975l : ISenguptal . 1 1 995b . Both cases can be regarded as 
integrable systems with the electromagnetic field acting as a perturbation. 

The above-mentioned lobes are defined by the figures of the effective potential 
in the poloidal plane. These were previously studied in the cont ext of charge sep - 
aration that is expected to occur in pulsar magnetospheres (e.g. Neukirchl . 1993 ). 
Here, we address whether the trajectories within these lobes are regular (i.e., whether 
the system is integrable), or if t hey instead exhibit a c haotic behavior. A related 
problem was studied recently by lTakahashi & Koyamal d2009) in an attempt to find 
a connection between chaoticness of the motion and the spin of a rotating black 
hole residing in the center. These authors suggest that chaotic behavior occurs for 
certain values of the black hole spin, while for others the system is indeed regular. 

The idea of investigating the connection between the spin of a black hole and 
chaoticness of motion of matter near its horizon is very interesting for the follow- 
ing reason. Because of high degr ee of symme try of the background spacetime, 
the unperturbed motion is regular ( Carter . 19681 ): no chaos is present. The electro- 
magnetic perturbation may trigger the chaos, however, its effect can be expected 
to diminish very near the horizon, where strong gravity of the black hole should 
prevail. This is also the region where the spin effects are most prominent. Further 
out various other influences become important due to distant matter and the tur- 
bulence in accreted material. Therefore the connection between the spin and the 
motion chaoticness is best applicable in the immediate vicinity of the black hole. 
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Figure 2.21: In the left panel we present a poloidal section of the selected isocon- 
tours of the effective potential Veff (^5 6). eq- (I2.73L for a charged particle (qQ = 2, 
L = 5 M) on the Kerr background (a = 0.5 M). We assume the presence of 
Wald uniform magnetic field (qBo = 2M~^). The off-equatorial potential lobes 
are present, allowing stable motion. Two exemplary trajectories of test particles 
are shown - in the left lobe a chaotic orbit of energy E = 1.796, while in the right 
lobe the regular, purely off-equatorial trajectory of ^ = 1.78. Both particles were 
launched at r(0) =3.11, 9(0) = 7i/4 with w^(0) = and their trajectories inter- 
weave with each other. We plot the poloidal (r, 9) projection of the trajectory; what 
appears as a lobe in the poloidal plane is an axially symmetric 3-dimensional ro- 
tational structure. The latter is illustrated in the right panel where the case of the 
off-equatorial regular trajectory is shown. 



i.e. within the inner parts of corona. 

The recurrence analysis ( Marwan et al. . 2007 1) provides us with a powerful tool 
for the investigation of complex dynamical systems. The method examines the 
recurrences of the system to the vicinity of previously reached phase space points. 
It has been typically adopted to study the experimental data, where often only some 
(if not just one) of the phase space variables are known from the measurements. 
Takens' embedding theorems ( Takens. . J981 ) are then used to reconstruct the phase 
space portrait of such a system. In our study we are equipped with the full phase 
space trajectory from the numerical integration of the equations of motion, so that 
we can use the recurrence analysis directly. 

It appears that the method of Recurrence Plots has not been employed in the 
context of relativistic astrophysical systems yet. To this end, one needs a consistent 
definition of the neighborhood of a point in the phase space in a curved spacetime. 
Below, we discuss the phase space distance and suggest a form of the distance 
norm suitable in such circumstances. 
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Figure 2.22: The overview of possible topologies of the off-equatorial potential 
structure above the event horizon (thick line in plots) of Kerr black hole endowed 
with the Wald test field. 



2.7 Equations of motion and the effective potential 

The phase space trajectories of integrable systems are regular, meaning that they 
are bound to the surface of an ^-dimensional torus, where n is the number of de- 
grees of freedom. The torus is determined uniquely by n constants of motion that 
are present in such a system. Its behavior can be explored by Poincare surfaces 
of section, which are defined by intersections of the ph ase space trajectory with a 
2-dimensional plane dLichtenberg & Liebermanl . Il992h . On the other hand, non- 
integrable chaotic systems generally have fewer integrals than the number of de- 
grees of freedom. In general, both the regular and the chaotic orbits may coexist in 
the phase space of a single system. 

Chaotic orbits are ergodic on the given hypersurface. Its dimension is now 
larger than n, and the section points thus fill areas in the plot of the Poincare 
surface. However, depending on the initial conditions, regular orbits can also 
appear in non-integrable systems. Such orbits maintain the value of some ad- 
ditional constant of motion, although it is not generally possible to write this 
constant in an explicit form. In the context of motion around black holes per- 
turbed by (weak) external sou r ces, various aspects of chaos were studied e.g. by 
Karas & Vokrouhhckvl (ll992|} j_lNakamura & Ishizukal d 1993b : IPodolskv & Veselv 
(ll998|), and very recently by[Sraierak & Sukova (2010). 

A standard approach to an integrable system with a non-integrable perturbation 
assumes complete control over the strength of the perturbation (i.e., the perturba- 
tion can be set to be arbitrarily weak). If this were the case, we could first switch 
the perturbation completely off, analyze the orbits, and then observe the impact 
of gradually increasing the perturbation strength upon these orbits. However, the 
class of off-equatorial bound orbits only exists when the electromagnetic term is 
strong enough to balance the gravitational attraction of the central body. Then the 
(sufficiently strong) perturbation is by itself the cause of the new kind of the regular 
motion that happens outside the equatorial plane. 
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Figure 2.23: Regular motion in the equatorial potential lobe in the fully integrable 
system of charged test particle (E = 0.99, L = 5M, ^ = 10^, r(0) = 32.02M, 9(0) = 
1.54) in the pure Kerr-Newman spacetime (Q = 3 x 10~^, a = 0.5 M) endowed 
with the fourth Carter constant of motion L . Long diagonals parallel to the LOI 
are general hallmark of regularity in the RPs. 



Having this delicacy on mind, we shall use the usual Hamiltonian formalism 
to express equations o f motion governing the trajectories. We first construct the 



super-Hamiltonian 9{ dMisner et al.l . ll973h . 



^ = y^'^^M - q\){^v - qAy), (2.69) 

where m and q are the rest mass and charge of the test particle, %^ is the generalized 
(canonical) momentum, g^^ is the metric tensor, and denotes the vector potential 
of the electromagnetic field. The latter is related to the electromagnetic tensor 
F^v by F^v =^v,;u ~^A^,v- Unless otherwise stated, we will use geometrical units, 
G = c = 1 (see Appendix lAb. 

The Hamiltonian equations are given as 

d^^ diz, 

where ?i = x/m is the affine parameter, x denotes the proper time, and is the 
standard kinematical four-momentum for which the first equation reads — Ti^ — 

In the case of stationary and axially- symmetric systems, we identify two con- 
stants of motion, namely, the energy E and angular momentum L. From the second 
Hamiltonian equation (12.701) we obtain 



T^t = A + q^t = -E 



(2.71) 
(2.72) 



CHAPTER 2. MOTION IN BLACK HOLE MAGNETOSPHERES 54 




r[M] time 



Figure 2.24: Regular off-equatorial motion of a charged test particle (E = 1.77, 
L = 5M, r(0) = 2.9M, 9(0) = 0.856 and u'{0) = 0) on the Kerr background (a = 
0.5 M) enriched with the Wald test field (qBo = 2M"\ qQ = 2). The diagonal 
structures typical for trajectories in integrable systems are preserved, though the 
pattern is more complicated. 



The trajectory is specified by the integrals of motion E, L, and the initial values 
r(0), 9(0) and w^(0). The initial w^(0) can be calculated from the normalization 
condition, g^^u^Uy = — 1 (we always choose the non-negative root). 

The effective potential can be derived in the following form: 

Veff= , (2.73) 

where 

a = (2.74) 

(3 = 2[g'^{L-qA^)-g^Atl (2.75) 

Y = -g'^'^^L-qA^f-g^'q^A^ + lg'^AtiL-qA^)-!, (2.76) 

and where we introduce specific quantities L=^,E = ^ and the specific charge 
q=^' Local minima of Veff (^5 6) reflect the location of stable orbits of test parti- 
cles. Off-equatorial potential minima were identif i ed and various types of potential 
lobes were discussed elsewhere (see Kovaf et al. . 2010 ^. We can express the ef- 



fective potential (12.731) also as a function of r and u^, and use it to determine the 
boundaries of allowed regions in Poincare surfaces of section for a given value of 
9. 



^We have employed the metho d of effective potential to study the stability of the m otion, however, 
we note that the force formalism (lAbramowicz et allll995l : [Kovar &Stuchli1<ll2007h can serve as a 
very efficient alternative tool. In particu lar, the off-equator ial motion of charged particles can be 
examined via the procedure described in Kovaf et This allows us to localize minima of 

the effective potential around which the stable orbits occur. 



CHAPTER 2. MOTION IN BLACK HOLE MAGNETOSPHERES 55 




Figure 2.25: A transitional state between the regular and chaotic regimes of mo- 
tion of a highly charged test particle which only differs from the previous case 
by increasing the energy to £ = 1.796. The diagonal lines in the RP are partially 
disrupted, indicating the onset of chaos. 

We employ the Kerr metric in standard Boyer-Lindquist coordinates r, 9, (p 
given by eq. (12.11) . It is sufficient to consider positive values of the spin parame- 
ter a without loss of generality (the cases of prograde and retrograde motion are 
distinguished by the sign of the particle charge and the orientation of the magnetic 
field). By setting a = the metric (12.11) reduces to the static one describing the 
Schwarzschild spacetime. 

At this point, a note is worth on the adopted computational scheme which 
we have employed to study the trajectories and to detect the chaotical behavior. 
In order to reach reliable results, we coded several approaches and we checked 
their stability and precision. We employ the multi-step Adams-Bashforth-Moulton 
solver to determine the phase- space trajectory by numerical integration of eqs. 
(12.701) . In some cases, when a higher precision is demanded, we use the 7-8th 
order Dormand-Prince method that belongs to the family of explicit Runge-Kutta 
solvers with adaptive stepsize. This method improves the accuracy significantly, as 
can be verified by checking the conservation of the integrals of motion along the 
trajectory. However, the improved accuracy comes at the expense of computational 
time, as the adopted Dormand-Prince scheme is more computationally demanding 
than the Adams-Bashforth-Moulton solver. 

Furthermore, it is well-know that, when dealing with Hamiltonian systems, the 
most appropriate solv ers are those w hich respect the symplectic nature of Hamil- 
tonian dynamics (e.g. l YoshidaL [l993h . We therefore employed also the implicit 
Gauss-Legendre Runge-Kutta (GLRK) method, which is a symplectic scheme. In- 
deed, we confirm that this code provides the most reliable results, especially in 
the case of long-term integration. The difference in the accuracy between GLRK 
and non- symplectic solvers reaches several orders of magnitude and it is generally 




Figure 2.26: The chaotic motion of a highly charged test particle which only differs 
from the previous case by increasing the energy to ^ = 1.7975. The diagonal lines 
in the RP are now disrupted and complex large-scale structures appear which are a 
characteristic indication of deterministic chaos. 



more apparent in the case of chaotic trajectories, as expected. However, the cost 
in terms of the computational time is also non-negligible, and so we only use the 
GLRK method to achieve very accurate long-time determination of the trajectory 
in several exemplary runs. Detailed comparison of the integrators is presented in 
Appendix El 



2.8 Recurrence analysis 

The Kerr metric is well-known and the analysis of test part icle motion in this space 



time was carried out in many papers (IMisner et al.Lll973h . Among important fea 



tures of the Kerr metric is the fact that the particle trajectories are integrable, and 
so the chaos can set in only when perturbations of the background gravitational 
field are introduced or additional electromagnetic interaction with fields of exter- 
nal sources are allowed. This is also where recurrence analysis can be helpful. 

Methods of phase space recurrences have been successfully applied to a wide 
range of various empirical data, not only in physics but also related to physiology, 
geology, finances and other fields. Recurrence plots are especially suitable for the 
investigation of rather short and nonstationary data. On the other hand, the method 
of recurrence analysis has not yet been widely applied to study the dynamical prop- 
erties of motion in relativistic systems. We thus briefly summarize this approach 
for our context. 

Besides more traditional methods of the numerical analysis of dynamical sys- 
tems, such as a visual s urvey of Poincare surfaces of section or the evaluation 



of the Lyapunov spectra iSkokosl . E010 ). the recurrence analysis is a rather novel 



technique, based on the analysis of recurrences of the system into the vicinity of 
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Figure 2.27: The test particle (L = 6M, qBo = and qQ = 1) is launched from 
the locus of the off-equatorial potential minima r(0) = 3.68 M, 9(0) = 1.18 with 
w^(0) = and various values of the energy E. In the left panel we set £ = 1.58 
and we observe ordered off-equatorial motion. For the energy of £ = 1.65 cross- 
equatorial regular motion is observed (middle panel). The trajectory occupies only 
a part of allowed potential lobe, regardless the length of the integration period. Fi- 
nally in the right panel with ^ = 1 .75 we observe irregular motion whose trajectory 
would ergodically fill whole allowed region after the sufficiently long integration 
time. We show that the motion is chaotic in this case. Spin of the black hole 
is a = 0.9 M and its event horizon is depicted by the bold line. Topology of the 
potential lobes corresponds to the type Id of our classification (see fig. 12.221) . 



its previous states. 

Recurrence Plots (RP) are int roduced as a tool of v isualizing the recurrences 



of a trajectory in the phase space (lEckmann et al.Lll987h . The method is based on 



examination of the binary values that are constructed from the trajectory x{t). Re- 
sults of the orbit analysis can be quantified statistically in terms of the Recurrence 
Quantification Analysis (RQA)|1 

The RP construction is straightforward regardless of the dimension of the phase 
space. We only need to evaluate the binary values of the recurrence matrix Rij, 
which can be formally expressed as follows: 

Rij{E) = 0(e- ||x(/) -x(j)ll), ij = 1, ...,A^, (2.77) 

where 8 is a pre-defined threshold parameter, the Heaviside step function, and 
N specifies the sampling frequency. The sampling frequency is applied to the time 
segment of the trajectory x{t) under examination. There is, however, no unique pre- 
scription for the appropriate definition of the phase space norm 1 1 . 1 1 in eq. (12.771) . 
We can consider a purely abstract vector space and apply one of the elementary 
norms L\ (Euclidean norm) or L°° (maximum norm). Some aspects of the ap- 
propriate choice of the norm are deferred to sec. 12.8.11 



^We use the CRP ToolBox (iMarwan et all l200l p. 321) in Matlab (R2009b) to construct RPs 
and to evaluate RQA measures. 
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Figure 2.28: Motion of the charged test particle (L = 6M, qQ = 1, qBo = M ^ 
r(0) = 3.68M, w"(0) = and 9(0) = Gsection = 1.18) on the Kerr background (a = 
0.9 M) with the Wald's test field is discussed. For E = 1.578 we observe regular 
off-equatorial motion (upper panels). Increasing the energy level to ^ = 1.65 we 
obtain cross-equatorial regular trajectory (middle panels). Rising the energy to 
^ = 1.75 we get a chaotic cross-equatorial orbit. In the Poincare surface of section 
we distinguish >0 (black point) from < (grey point). Analyzed trajectories 
are of the type Id (fig. 12.221) and correspond to those of fig. 12.271 
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Figure 2.29: Color map of the mutual phase space distance seperating the given pair 
of points on the trajectory. This example concerns a charged particle trajectory near 
the Kerr black hole in the asymptotically uniform magnetic field. Left: the case of 
regular motion with the energy of ^ = 1.77. Right: the case of chaotic motion with 
E = 1.7975. In the latter case more complex structure appear in the recurrence 
plot. The common parameters of both panels are L 5M, a 0.5 M, qBo 2M~^ 
and qQ = 2 with the initial condition r(0) = 2.9M, 9(0) = 0.856 and w"(0) = 0. 



Finally, we need to speci fy the value of the thr eshold parameter 8. To this end 
we follow the suggestion of Marwan et al.l (120071 . sec. 3.2.2) and relate 8 to the 
standard mean deviation, a, of the given data set. Setting 8 = is advantageous 
because the proportionality constant k, once adjusted to obtain a properly filled Re- 
currence Plot, remains valid (with only minor adjustments) for all data sets of other 
trajectories. We therefore normalize the time series of each coordinate separately 
to zero mean and o = 1. 

The binary valued matrix R^^ represents the RP which we get by assigning a 
black dot where Rij = 1 and leaving a white dot where Rij = 0. Both axes represent 
a time segments over which the data set (the phase space vector) is being examined. 
RP is thus symmetric; the main diagonal is always occupied by the line of identity 
(LOI). 

Recurrence Plots co ntain wealth of information about the dynamics of the sys- 
tem (ThieletaL, 2004ab . Different pieces of k nowledge are encoded in large-scale 
and short-scale patterns (IMarwan et al.L l2007l . sec. 3.2.3). To decide if a particular 
trajectory is a regular or a chaotic one, the determining factor is the presence of 
diagonal structures in the RP. Diagonal lines in the RP reflect the time segments 
of phase space trajectory during which the system evolution proceeds in a regu- 
lar way. It captures the epoch when the trajectory proceeds almost parallel to its 
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Figure 2.30: Graphs of different RQA measures based on the diagonal lines in the 
RP as a function of specific energy E. In each panel, 400 trajectories were ana- 
lyzed in a given energetic range. All measures exhibit the evident change of their 
behavior aiE ^ 1.7954 which we interpret as an onset of chaos. Other parameters 
remain fixed at following values: L = 5M, qBo = 2, r(0) = 2.9M, 9(0) = 0.856, 
w"(0) = 0, a = 0.5M and qQ = 2. 



previous segment, i.e. within the 8-tube around that segment. Hence, integrable 
systems exhibit themselves by diagonally oriented structures in their RR On the 
other hand, if the motion is chaotic the diagonal lines disappear and the diagonal 
features become shorter, as the trajectories tend to diverge quickly. As a result, 
more complicated structures appear in the RR 

We stress th at the interpretation o f the RP is primarily intuitive. We refer to 
the review paper iMarwan et al.l (120071 . sec. 3.2.3) where the patterns appearing in 
the RP and their relation to the current dynamic regime are analyzed in detail. 
Although basic conclusions may be inferred in general (e.g. distinction between 
regular versus chaotic regime) the fine structure of the RP depends heavily on the 
properties of a given dynamic system. In order to gain more insight into the way in 
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Figure 2.31: Shannon entropy of probability distribution of diagonal lines lengths 
ENTR and three RQA measures based on the vertical lines of RPs as a function of 
specific energy E (details in the text). To some surprise, the vertical measures also 
react dramatically to the onset of chaos at ^ = 1.7954. 

which various dynamical regimes manifest themselves in our system we typically 
present RPs accompanied by corresponding Poincare surface of section throughout 
this text. 

Visual behavior of RP and its complexity is quantitatively reflected in RQA. 
The RQA evaluates statistical characteristics of the recurrence matrix R/^. First of 
all, we define the recurrence rate (RR) as a density of points in RP, 

'^R(^) ^Wit ^'-^■(^)- (2.78) 

Now we can turn our attention to diagonal segments in RP. Their length draws 
distinction between regularity and chaos. The histogram P(8, /) records the number 
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of diagonal lines of length /. It is formally given as follows: 

P(8,/) = £ (1 -R,_ij_i(8))(l -R,+/j+,(8)) nR,+,j+,(8). (2.79) 

i,j=\ k=0 

This histogram defines the determinism factor (DET), defined as a fraction of recur- 
rence points, which form the diagonal lines of length at least Imin to all recurrence 
points, 

DET=-= — ^ — -. (2.80) 

The average length of diagonal lines L (where only lines of length at least Imm 
count) is 

Lf-f. lP{z,l) 
L=-= — ^ — -, (2.81) 
L^f. P{e,l) 

and the corresponding divergence (DIV) is defined as inverse of the length of the 
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Figure 2.33: Shannon entropy of probability distribution of diagonal lines lengths 
ENTR and three RQA measures based on the vertical lines of RPs as a function of 
specific energy. Sudden change of the behaviour at E ^ 1.685 is apparent for all 
the quantities. 



longest diagonal line Ln 



DIV 



(2.82) 



DIV is in its very nature closely rela ted to the divergent fea tures of the phase space 
trajectory, and so it was originally (lEckmann et al.L Il987h claimed to be directly 
related to the largest positive Lyapunov characteristic exponent ?imax- On the other 
hand, theoretical considerations justify the use of DIV as an estimator only for the 



lower limit of the sum of the positive Lyapunov exponents dMarwan et al.L 120071 . 
sec. 3.6). Neverthel ess, a strong correlation between DIV and ?Lniax arises in nu- 
merical experiments ( TruUa et aP . [7996 ). 



The quantification measure ENTR is defined as the Shannon entropy of the 
probability p(8, /) = P(8, /) /Ni of finding a diagonal line of length / in the Recur- 
rence Plot, 



ENTR=- In ;?(£,/), 



(2.83) 



where Ni is a total number of diagonal lines: Ni{e) = lli>i^i^ 0- 

Analogous statistics may be performed for vertical as well as the horizontal 
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Figure 2.34: Comparison of purely off-equatorial trajectories in spacetimes differ- 
ing by the spin parameter a (left panels: a = 0.3M; middle: a = 0.6M; right: a^M) 
which is linearly linked to the energy E (left panels: E — 1.56; middle: E — 1.9, 
and E = 2.35 in the right panels). Other parameters remain fixed: L = 5M, 
6(0) = Osection = 0.856, qB^ = 2M~^ and qQ = 2. RPs are taken for trajectories 
with r(0) = 2.9 M and w^(0) = 0. Although the structures in the Recurrence Plots 
clearly differ from each other, all of them represent diagonally oriented patterns 
that are characteristic of regular motion. The regularity of the motion is confirmed 
by surfaces of section in the bottom panels, where several trajectories (for each 
value of spin and energy) are presented. Different colours are used to distinguish 
the orbits originating from different initial conditions. 



segments (RP is symmetric with respect to the main diagonal). These segments are 
generally connected with periods in which the system evolves during its laminar 
state. To this end, the histogram P(£,v) records the number of vertical lines of 
length V and it can be constructed as follows: 

P(8,v) = f (1 -R,-Xe))(l -Ru+v(e))nRu+^(e). (2.84) 

In analogy with the diagonal statistics histogram, P(£,v) is used to define the 
vertical RQA measures. Laminarity (LAM) is defined as a fraction of recurrence 
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points that form vertical lines of length at least Vmin to all recurrence points, 

LAM = t^^^ — 1^. (2.85) 

The trapping time (TT) is an average length of vertical lines, 

TT = '^,7^°"" — . (2.86) 

Finally, the length of the longest vertical line (Vmax) can also be of some interest. 

RQA measures the crucial dependence of RP on the value of the threshold 
parameter, 8, which must be adjusted appropriately to a given data set. This lack 
of invariance is a drawback of both RPs and RQA. Nevertheless, it was shown 



(IThiel et al.L |2004b) that stable estimates of various dynamical invariants, such as 



the second order Renyi entropy and the correlation dimension, can be inferred if 8 
is kept within a reasonable range. Since we shall use the standard RQA measures 
to compare the dynamics between test particles with different initial conditions, 
we have to eliminate the numerical effect of variances in the range of coordinate 
values spanned by these trajectories. We achieve this by fixing the value of 8. 

After the brief set of preliminaries we are now prepared to proceed to the in- 
tended application of the recurrence analysis in the next section. 

2.8.1 Phase space recurrence in the curved spacetimes 

When analyzing the dynamics in the general relativistic context the fundamental 
question arises whether the distinction between chaotic and regular motion is co- 



ordinate dependent or not. To this end iMotterl (l2003h infers the transformation 



law for Lyapunov exponents. He concludes that although the Lyapunov exponents 
themselves are not invariant they transform in such a way that positive Lyapunov 
exponents remain positive and vice versa. In other words the distinction between 
regular and chaotic dynamics may be drawn invariantly. 

In order to give the notion of recurrence a rigorous and, at the same ti me, an 
intuitive sense, we can employ the 3 + 1 formalism (Tho me & MacDonaldl . U982) 
that is based on an appropriate selection of a family of spacetime-filling three- 
dimensional spacelike hypersurfaces (foliation) of constant time t. The timelike 
curves orthogonal (in a spacetime sense) to the hypersurfaces may be regarded as 
the world-lines of a family of fiducial observers (FIDO) who naturally parameterize 
their world line by proper time x (whose rate of change generally differs from that 
of t). FIDO identify each spatial hypersurface along his world line as a slice of 
simultaneity. The geometry of this spacetime slice is given by 3 -metric Y/j^ 

lij = gij + mj, (2.87) 

where w/ stands for the spatial part of FIDO's four- velocity and gtj for the spatial 
part of the spacetime metric. Considering also the time coordinate, y^v can be 
regarded as a projector onto the three-dimensional spatial hypersurface. 
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Figure 2.35: Comparison of trajectories of particles launched from the equatorial 
plane with different spin values. Also in this case we have to link linearly the value 
of spin a with E in order to maintain the existence of the potential lobe. In left 
panels we set a 0.5M, E 1.795, in middle panels a 0.6M, ^ = 1.92 and 
in right panels a — M, E — 2.42. For all three cases we show surfaces of section 
of several trajectories differing in initial values r(0) and w^(0). Recurrence Plots 
are taken for trajectories with r(0) = 2.15M, w^(0) = 0. Other parameters remain 
fixed: L = 5M, qB^ = 2M-\ 9(0) = Gsection = f , ^2 = 2. 



Our selection of FIDOs will be restricted here to those moving along the Killing 
direction (thus seeing an unchanging spacetime geometry in their neigbourhood). 
We shall also demand that all the hypersurfaces of simultaneity which they con- 
stitute have an identical spatial geometry. Finally, we require that FIDO's proper 
time T becomes asymptotically identical to the coordinate time t. In the case of 
Kerr spacetime (as an example of stationary and axisymmetric spacetime) these 
restrictions lead to the identi fication of the FIDOs with th e zero angular momen- 
tum observers (ZAMOs), see lThome & MacDonaldl (|l982l ). Orthonormal tetrad of 
ZAMO is given by eqs. ([244l) - ([247]) . 

Projecting an arbitrary four- vector onto ZAMO's hyper surf ace of simul- 
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2.5 




Figure 2.36: Profiles of the effective potential Veff for L = 5M, a = 0.5 M, qQ = 
1.03, qBo = 0.1M~^, taken in the equatorial plane and in the asymptotic region 
z = rcosO oo (as indicated with the corresponding curves). The valley crosses 
the equatorial plane, almost unaffected in its bottom parts, although the behavior of 
the potential near the symmetry axis is quite different in the equatorial plane where 
it approaches the horizon of the black hole (vertical dotted line at r = r+(a)). The 
horizontal dotted line at unity measures the rest energy of the particle m. 



taneity directly results in the following 3-dimensional quantities, 

= ifC^=(C„Ce,gcpcpC^); 

here, Ycpcp = gcpcp (W(p = for ZAMO) and, consequently, YP^g(p(p = 1. 
ZAMO tetrad components ^^C^^^ are given as 



(2.88) 
(2.89) 



(2.90) 
(2.91) 



The hypersurface components of the phase space constituents x' and Ji,-, as mea- 
sured by ZAMO, are then: 



1 



3D 



The spatial 3-metric is 



(2.92) 
(2.93) 

(2.94) 
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where x^^^ represents the spatial distance from the origin of the tetrad, i.e. ZAMO's 
current location. ZAMO is not an inertial observer, which generally causes the 
first order corrections 0(|x^^^|) to the Minkowskian metric g(i)(j) = ^ (/)(;)• 
these do not enter the spatial part of the metric. Thus the 3-metric within the 
spatial hypersurface is a Euclidean one, with the deviations of second order in the 
distance from the spatial origin on ZAMO's world-line. 




time time time 



Figure 2.37: An exemplary trajectory (E = 1.058, L = 5M) is launched from 
the equatorial plane 9(0) = | with w^(0) = 0. Parameters of the background are 
a = 0.5M, qBQ = 0.1M~\ qQ = 1.03. In the upper left panel we observe that set- 
ting r(0) = 8.4M results in oscillations around the equatorial plane while launch- 
ing it at r(0) = 8.7M makes it escape. In the upper middle panel we examine 
the trajectory of the escaping particle in terms of the rescaled radial coordinate 
r* = ^-y^. In the case of oscillating trajectories two distinct modes of motion are 
observed (upper right panel). The first particle (r(0) = 11. 5M) shows a complex 
"ribbon-like" trajectory; the other one (r(0) = 8.4M) fills uniformly the given por- 
tion of the potential valley. The Recurrence Plots are also shown (bottom panels). 
We observe a highly ordered regular pattern for the particle with r(0) = 8.4M (left 
panel), a more complicated diagonal pattern of the ribbon-like trajectory (launched 
at r(0) = 11. 5M, middle panel), and a disrupted diagonal pattern of the transitional 
trajectory (r(0) = II AM, right panel). 

We will use ZAMO's metric at distances up to the value of the threshold pa- 
rameter 8. The Euclidean metric according to eq. (12.941) will be therefore justified 

2 

if ^ <^ I, where P stands for a constant (for a given ZAMO). P has the dimen- 
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sion of length and characterizes the curvature of the hypersurface. We suggest 
setting P = K~^/^, where K = T^^^^^/^^y^^^ represents the Kretschmann scalar eval- 
uated from the Riemann curvature tensor. In the case of the Kerr black hole th e 
Kretschmann scalar may be expressed in a surprisingly simple form (|Henry|,|2000|). 

2 

While constructing the Recurrence Plots, we check whether the condition ^ 
remains satisfied. 

The above-mentioned adoption of preferred observers is needed in order to 
maintain an operational criterion of chaos and be able to formulate an explicit 
form of the equations for RQA measures in a curved spacetime. Here the notion of 
the phase space distance plays a role. In Kerr metric (or another axially symmetric 
stationary spacetime), Fiducial Observers (a.k.a. FIDOs) represent a natural selec- 
tion of preferred observers. Obviously, this option is not unique, and so a detailed 
appearance of the recurrence plots is also ambiguous to certain extent. But not so 
the main conclusions that we infer regarding the chaoticness versus regularity of 
the system behavior, because this distinction can be eventually traced down to the 
exponential versus polynomial growth of the separation with the particle proper 
time along neighboring trajectories. 

We can deduce the kind of transformation between different families of ob- 
servers that could affect our conclusions: these are transformations involving ex- 
ponential dependencies on observer's phase-space position. For example transfor- 
mation to accelerated frames and spacetime points in the vicinity of singularities 
may need a special consideration, as well as the investigation of highly dynami- 
cal spacetimes that are lacking symmetries. On the other hand, selecting LNRF to 
define ZAMOs in (weakly perturbed) Kerr metric outside the black hole horizon 
appears to be a well-substantiated choice. 

Similar arguments for the adoption of preferred observe rs on the basis of space- 
time s ymmetries have been elaborated in greater detail by iKaras & Vokrouhlicky 
(I1992I ) in the context of Ernst's magnetized black hole, which is another partic- 
ularly simple (static) exact solution of Einstein-Maxwell equations exhibiting the 
onset of chaos as the magnetic field strength is increased. 



2.9 Kerr black hole in uniform magnetic field 



Large-scale magnetic fields are known to be present in cosmic conditions. They 
can exist around black holes, which do not support their own magnetic field but 
may be embedded in fields of distant sources. In the case of neutron stars, dipole- 
type magnetic fields of very high strength often arise. We concentrate on black 
holes in this section and defer the case of a magneti c star to sec . 12.101 

By employing the uniform test field solution dWaldl . 1 19741 ) we incorporate a 
weak large-scale magnetic field near a rotating black hole. The vector potential 
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Figure 2.38: The Poincare surface of section of several trajectories (E = 1.058, 
L 5M, qBo 0.lM-\ qQ 1.03, a 0.5M, 9 = f ) launched from the equatorial 
plane with various values of r(0) and w^(0) = 0. Grey colour indicates the escape 
corridor which lets the particles escape from the equatorial plane (see details in the 
text). 



can be expressed in terms of Kerr metric coefficients (12.11) as follows. 

At = ^Bo{gtc^ + 2agtt)-^Qgtt-\Q. (2.95) 
A(p = \Bo (g(p(p + 2agti^) - \Qgti^, (2.96) 

where is magnetic intensity and Q stands for the test charge on the background 
of Kerr metric. The terms containing Q can be identified with the components 
of the vector potential of Kerr-Newman solution (although the test charge does 
not enter the nietric it self). An example of an integrated trajectory is shown in 
fig. [2211 IWaldl (119741 ) has shown that the black hole selectively accretes charges 



from its vicinity, until it becomes itself charged to the equilibrium value 

(2w = 25oa. (2.97) 

We remark that the particle charge q appears always as a product with Q or 
in the formula (12.731) for the effective potential, as well as in equations of motion 
(12.701) . Therefore, the simultaneous alteration of q, Q and 5o values, preserving 
the products qQ and qB^, does not affect the particle dynamics. If we further 
assume that Q = = 25o<2 is maintained, we only need to specify the value of 
qB^ to uniquely determine a particular trajectory. However, since we do not restrict 
ourselves to the case Q = we decide to always explicitly specify the values of 
qQ and qB^. 
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Figure 2.39: Selected types of the effective potential behavior in the vicinity of 
off-equatorial halo orbits above the surface of magnetic star with rotating dipole 
magnetic field. The inner bold line signifies the surface of the star at r = 4M. The 
outer line is the light surface. 



2.9.1 Motion within the potential lobes 

Our previous analysis (Kovaf et al.l . l20ld ) concludes that the off-equatorial bound 
orbits are allowed only for test particles obeying simultaneously the two conditions, 
sgn(aL) = 1 and sgn(^) = sgn(a5o). Four distinct types were found (see fig. 12.221) 
and we examined the dynamics of test particles in all of these types. The results 
for different types are comparable, and so we present here only the analysis of two 
of them, namely the types la and Id. While raising the energy level from the local 
minima of the symmetric halo orbits, we observe that the off-equatorial lobes grow 
and eventually merge with each other once the energy of the saddle point in the 
equatorial plane is reached. If we further increase the energy the behaviour between 
the two classes differs profoundly. For la we observe that the lobe eventually opens 
toward the horizon allowing the particle to fall onto the black hole. On the other 
hand in the Id setup the merged lobe breaks "toward infinity" allowing the particle 
to escape in the axial direction when the energy is raised sufficiently. 

The size of the lobes is controlled by the specific energy E. Employing the 
Poincare surfaces of section and the Recurrence Plots we investigate how the regime 
of the particle motion changes with E. This parameter appears as a suitable control 
parameter producing a sequence of bound trajectories while all other parameters 
(and the initial position) remain fixed. For the sake of comparison we first present 
the case of a fully integrable system of a charged test particle on the pure Kerr- 
Newman spacetime (fig. 12.231) . The motion occurs in the potential lobe around the 
local minim um in the equatorial plane; there ar e no halo orbits above the horizon 
in this case (iKovaf et al.Ll20ld : lde Felicel . ll979b . 

On the Kerr background with Wald test field we first analyze the sequence 
of trajectories of la type. Figure 12.241 shows regular motion occurring in the off- 
equatorial lobe. Increasing the energy level (while keeping all other parameters 
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fixed) above the value in the equatorial saddle point results in a transitional regime 
depicted in fig. 12.251 The onset of chaotic features does not occur as a direct con- 
sequence of the lobe merging. We rather observe that the orbit bound in merged 
(cross-equatorial) lobe remains regular until the particle notices the possibility of 
crossing the equatorial plane which happens when its energy is increased suffi- 
ciently. Once the motion becomes cross-equatorial, chaotic features appear. By in- 
creasing the energy even more we approach the critical value when the lobe opens 
and allows the particles to fall onto the horizon. For energies slightly below this 
limit we detect a fully chaotic regime of motion (fig. 12.261) . 

Similar sequence of trajectories is analyzed in the case of Id class of the topol- 
ogy of the potential wells. Three typical trajectories found in this case are shown: 
regular off-equatorial, regular cross-equatorial and chaotic cross-equatorial orbit. 
First we present their overview in fig. 12.271 which shows the main difference from 
the la series clearly: the presence of fully regular cross-equtorial orbits. In fig. l2.28l 
we compare Poincare surfaces of section and Recurrence Plots of these trajectories. 
Unlike previous series we distinguish >0 from < in the Poincare surface 
of section. 

Figure 12.291 shows an alternative representation of the recurrence plots, where 
the phase space distance separating the given pair of points of the trajectory is 
encoded by different colors. Again, by comparing the two panels one can clearly 
recognize how the diagonal structures disintegrate into scattered points as the chaos 
sets in. 

From the survey of the first series of Poincare surfaces of section and the Re- 
currence Plots (figs. l2.24H2.26l) we may only conclude that the transition from the 
regular to chaotic regime occurs somewhere close to the value E = 1.796. In order 
to localize this transition more precisely we evaluate RQA measures for 400 trajec- 
tories with the energy spread equidistantly over the interval E G (1.7948, 1.7968). 
In figs. (12.301) and (12.311) , we observe a sudden change of the behavior of the statis- 
tical measures at £ ^ 1 .7954, reflecting a dramatic change of the particle dynamics. 

Moreover, we know that the divergence DIV is related to the Lyapunov expo- 
nents, and in fig. 12.301 we observe that it suddenly rises at E ^ 1.7954, meaning 
that the trajectories become more divergent when this energy is reached. All of 
these indications combined lead to the conclusion that this energy level represents 
a critical value at which a transition from regular to chaotic regime occurs. 

Similarly we perform the RQA for the sequence of orbits of Id type which 
we analyzed qualitatively in fig. 12.281 Visual survey suggests that the suspicious 
interval in which the transition from regular to chaotic dynamics shall occur is 
E G (1.65, 1.75). We calculate 200 trajectories with energies equidistantly spread 
over the given range while other parameters of the system are fixed at values used 
in fig. 12.281 In figs. 12.321 and 12.331 we observe that all queried RQA measures 
exhibit sudden change in its behaviour ai E ^ 1.685. This is where the dynamic 
transition between the regimes occurs. Transition from the regular motion to the 
chaos is detected not only by diagonal RQA measures but also by the vertical ones. 

We conclude that the energy of the particle E acts as a governing factor deter- 
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mining the dynamic regime of motion. Our survey across various initial conditions 
has shown that motion in potential wells of the type la in a Wald test field is gen- 
erally regular. Chaos appears well after the merging point of the lobes and close to 
the critical breaking energy. We have verified that all queried RQA measures, i.e. 
not only those based on diagonal lines in RP, react to the transition from the regular 
to chaotic regime of motion, allowing us to localize precisely the transition. This 
was fully confirmed also in the other two types (classes lb and Ic of our typology 
in fig. [222]). 



2.9.2 The effect of spin on the chaoticness of motion 

Much attention has been recently focused towards the problem of detern iining the 
black hole spin from the proper ties of motion of surrounding matter (INarayanl . 
20051 : [Reynolds & NowakL l2003h . The astrophysical motivation to address these 
issues arises from the fact that cosmic black holes are fully described by three 
parameters - mass, electric charge and spin. While the methods of mass de- 



termination have be e n wid ely discussed (e.g. ICasaresl . l2007l : IVestergaardl . l2010l : 
Czemy & Nikolajuk . 2010b . the electric charge is considered to be negligible be- 
cause of rapid neutralization of black holes via selective accretion. However, de- 
termining the spin is a much more challenging task: the sp in is importan t , but its 
influence is apparent only very near the black hole horizon ( Murphy et al. . 2009b . 




Figure 2.40: Regular motion in the off-equatorial potential lobe at the energy level 
E = 0.8482. Parameters of the system are q^M = -5.71576 M^, L = 0.87643 M, 
Q = 0.011485 M~^. The left panel shows sections of several trajectories launched 
from 9(0) = Gsection = 1.0492. One of them (r(0) = 5.02 M and w"(0) = 0) is 
visualized in the Recurrence Plot in the right panel. The motion is regular (RP 
remains diagonal), however, the density of recurrence points clearly grows during 
the analyzed period. 



One can raise a question of whether the value of spin parameter a of the Kerr 
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black hole affects the dynamical regime of motion in the immediate neighborhood 
of the black hole. In other words, we ask if the spin parameter a triggers or dimin- 
ishes the chaoticness of the system. Answering this question is not straightforward 
because by altering the spin across an interval of values (a^ < M^) we inevitably 
have to change some other variables of the system; otherwise the different cases 
could not be directly compared. Moreover, the location and the very existence 
of the potential lobes is not automatically ensured over the whole range of spin 
because of strong Veff(^) dependence. 

We found that in order to keep the off-equatorial lobe at the initial position and 
the original size we need to increase the energy E, roughly proportionally to the 
increment of a. The effect of increasing a exhibits itself by lifting the hyperplane of 
effective potential. To compensate for this effect we have to elevate the ^-plane at 
which we cut the potential, so that we obtain (roughly) the original closed contour 
(the potential lobe) inside of which the motion is confined. It turns out that this can 
be achieved by linking both quantities linearly. 

First we compare the dynamics in the off-equatorial lobe in the range ^ E 
(0.3, 1) (for a < 0.3M the topology of the effective potential changes) to which 
we linearly relate the energy range E G (1.56,2.35) (whilst other parameters are 
kept fixed as follows: L = 5M, qBo = 2M-\ r(0) = 2.9M, 9(0) = 0.856, w"(0) = 
0, qQ = 2). By inspecting the Poincare surfaces of section and performing the 
recurrence analysis for a large number of trajectories across the given range of a 
and E (exemplary cases presented in fig. l2.34l) we come to the conclusion that there 
is no overall trend that could suggest that a is the unique driving agent affecting the 
regime of motion. All trajectories in our survey exhibit a regular behavior, which 
is also in agreement with the previous conclusion that the motion in off-equatorial 
potential lobes associated with the Wald test field is generally regular. 

We also examined the dynamics of test particles launched from the equatorial 
plane whose trajectories occupy the potential lobe extending symmetrically above 
and below the equatorial plane. A given lobe maintains its size for spin values ^ G 
(0.5, 1) and the related interval of energy E G (1.795,2.42). A survey across the 
given range of spin (energy) values reveals for this class of trajectories both chaotic 
and regular regimes. In fig. 12.351 we observe that for the lowest inspected spin, 
a = 0.5M (E = 1.795), the regular motion dominates, although islands of chaotic 
behavior are also present. Increasing the spin (energy) we observe that regular 
trajectories gradually diminish. For a = 0.6M (E = 1.92) some regular orbits still 
appear, but they are already dominated by chaotic trajectories. For higher spins 
the traces of regular motion further diminish. We present the extreme case a^M 
(E = 2.42) in fig. l2.35l to illustrate this apparent chaotic takeover. 

We conclude that for the class of orbits originating in the equatorial plane, spin 
a could possibly act as a destabilization factor which triggers the chaotic motion if 
enhanced sufficiently. However since the energy E is increased simultaneously it 
is not possible to attribute the observed dependence to the spin itself. On the other 
hand we have already seen in the above-given discussion (sec. 12.9.11) that energy E 
may itself act as a key factor determining the dynamic regime of motion. Thus we 
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Figure 2.41: For an energy value of ^ = 0.8485, both off-equatorial lobes merge 
via the equatorial plane. The upper-left panel shows Poincare sections of two tra- 
jectories (9(0) = Gsection = 1.0492, u'{0) = 0). A particle launched at r(0) = 4.8 M 
never crosses the equatorial plane and moves regularly. Setting r(0) = 5 M we 
observe a chaotic motion crossing the equatorial plane repeatedly. All particles 
launched with r(0), w^(0), corresponding to the inner parts of the potential curve, 
move in the same chaotic manner. The outskirts are occupied by regular trajec- 
tories. The upper-middle panel shows the transient trajectory (r(0) = 4.85 M), 
regular during the integration period of ?i = 10^. In the upper-right panel, the in- 
tegration time is prolonged to ?i = 3 x 10^. Here, the onset of chaos is connected 
with the first passage through the equatorial plane. The RP of the regular trajectory 
with r(0) = 4.8 M is presented in the bottom-left panel; the RP on the right belongs 
to the chaotic trajectory with r(0) = 5 M. 



suggest to attribute the observed triggering of the chaos to the increase of energy 
E rather then spin a. 

We remark that a similar problem conc erning the spin dependence o f the mo- 
tion chaoticness was addressed recently by Takahashi & Koyamal (2009). Authors 
of the quoted paper employ a dipole magnetic test field upon a Kerr background 
and perform a study of test particle trajectories, concluding that increasing the 
value of the spin parameter stabilizes the motion in a given setup. Unlike our case, 
the topology of the effective potential in their scenario allows the chosen potential 
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lobe to be maintained at a given location and roughly the same size (but not the 
same depth) even if a varies while E is kept constant. However increasing the spin 
allows the authors to set gradually lower and lower energies, for which more regu- 
lar trajectories are found. This is not at all surprising in perspective of our results, 
where the energy E proved to play a key role in determining the stability of motion. 
Although a direct comparison of presented surfaces of section differing only in a 
value may suggest the spin dependence, there is no clear and unambiguous corre- 
lation. We thus suggest attributing the observed dependence primarily to the level 
of energy, which acted as a motion destabilizer also in our setup. In fact, even if a 
real trend with the spin is present, it is hard to disentangle it from a simultaneous 
change of E. 

The question of the dependence of the dynamics upon the other parameters of 
the system (L, qBo, qQ) was also addressed during the analysis. It appeared that 
above mentioned difficulties accompanying the analysis of the spin dependence 
became even more serious in this case. Namely, neither it was possible to maintain 
the given potential lobe for a reasonable range of values of selected parameter nor 
we were able to fix this problem by binding this parameter in some simple manner 
to some other parameter (e.g. energy E). In other words, none of these parameters 
itself may be regarded as a trigger for chaos. 

2.9.3 Motion in potential valleys 

Besides off-equatorial lobes, the potential may form another remarkable structure 
- an endless potential valley of almost constant depth which runs parallely to the 
symmetry axis (fig. 12.361) . The poloidal orientation and asymptotical form of the 
valley are due to the fact that the Wald test field does not vanish at spatial infinity 
and it approaches the uniform magnetic field parallel to the symmetry axis. 

The existence of such a potential corridor suggests that test particles with a 
particular range of parameter values could escape from the equatorial plane to large 
distance. We observe that a test particle in the potential valley can keep oscillating 
around the equatorial plane or it may escape from this plane completely, depending 
on its initial position in the phase space (see the upper left panel of fig. 12.371) . We 
can examine the motion in the asymptotic region by rescaling the radial coordinate 
(upper middle panel of fig. 12.371) . To achieve this, we use r* = where r+ = 
M + \/M^ — d^ is the position of the outer horizon of the black hole. 

The normalization condition allows one to express as a function of other 
phase space variables and parameters of the system. Intuition suggests that, con- 
sidering particles launched from the equatorial plane (6(0) = f ), the initial value 
w^(0) is a governing parameter which decides whether the particle remains oscil- 
lating around 9 = | or leaves it once forever. We can draw the isolines of selected 
\^ values in the (r, w^)-plane which we use as a surface of section (Bsection = f ) 
for the inspection of the test particle dynamics. Comparing acquired isolines for 
various values of \^ with the empirically stated escape corridor of fig. [238] leads to 
the conclusion that they never coincide perfectly, although the correlation is quite 
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high. In other words there is no definite threshold value of w^(0) which would 
determine whether a selected combination of r(0), w^(0) (while other parameters 
are fixed) lets the particle launched at equatorial plane leave or oscillate around the 
plane. 

In fig. 12.381 we observe several qualitatively different types of possible particle 
dynamics. Next to the effective potential contour we find closed and well-defined 
curves which represent the regular motion. Going further inside we notice that 
fragmented curves are present. Below them we find progressively more and more 
blurred patterns, which again signifies the onset of chaos. The inner parts of the 
potential lobe are occupied by the "escape corridor" where the particles can stream 
freely from the equatorial plane, as seen in fig. 12.371 

The trajectories represented by the closed curves in the surface of section differ 
profoundly from the disconnected curve orbits. The difference can be seen also in 
the direct projection onto the poloidal plane (upper right panel of fig. 12.371) . While 
the trajectories of the first type gradually fill each particular compact region of 
given section, the latter forms bundles which curl through the projection plane 
resembling ribbons that bound regions which are never reached by the particle. 

By employing the Recurrence Plots (bottom panels of fig. 12.371) we confirm 
that the dynamics differs significantly in these distinct modes of motion. Not sur- 
prisingly we obtain a typical regular pattern in the case of trajectory which forms 
a closed sharp curve in the surface of section (bottom left panel of fig. 12.371) . 
A "ribbon-like" trajectory (fragmented curve in fig. 12.381) results in an ordered 
checkerboard pattern (bottom mi ddle panel), which is known to be typical for peri- 
odic and quasi-periodic systems (IMarwan et al.Ll2007h . Finally, in the bottom right 
panel of fig. 12.371 we observe that a blurred curve trajectory exhibits slight chaotic 
behavior in its RP. The diagonal structures are partially disrupted and we notice 
that diagonal lines become bent as they approach the line of identity. 

The cl ass of escaping trajectories in the Kerr background was recently dis- 
cussed by IPretil (120 id ). The author suggests that a Wald electromagnetic field 
employed in this setup could serve as a charge separation mechanism for astro- 
physical black holes since the sign of the particle charge may determine whether 
a given particle escapes from the equatorial plane or becomes trapped in cross- 
equatorial confinement (or falls into the horizon). 



2.9.4 Rotation number and fragmented tori 

In the integrable system the trajectories in the phase space reside on the surface 
of the tori characterized by the values of integrals of motion which determine the 
characteristic frequencies of the orbit. Fundamental frequencies in our axisym- 
metric system are those of radial and latitudinal motion ((p is cyclic coordinate). 
Once the system is slightly perturbed (by the magnetic field in our model) the 
tori characterized by the irrational ratio of frequencies co^/coe survive which is 
assured by the Kolmogorov-Arnold-Moser (KAM) theorem. On the other hand 
the Poincare-Birkhoff theorem tells us that the resonant tori with the rational fre- 
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Figure 2.42: Fixed point uq is localised. Two succeeding section points ui and 
of a given trajectory are marked. The angle 9/ between the radius vectors r/ 
and r/+i is evaluated. Rotation number v is defined as V = limA^^oo 2^ Y4L1 6/- In 
practise we compute a mean value of the finite number of angle values 0/ instead 
of the limit. 



quency ratio will disintegra t e into the chain of islands when we perturb the system 
(ILichtenberg & Liebermanl . Il992h . Trajectories belonging to the given chain are 
characterized by a single frequency ratio. Since the chains have nonzero radial 
width (whose actual value depends on the 'degree of prominence' of the resonance 
as well as the values of other parameters of the system) we should be able to detect 
them in a plot of the ratio 00^/ CO0 as a function of the initial radial coordinate of the 
orbit r(0). They should appear as the periods of constancy in such a plot. 

Rotation number is yet anoth er index which can be used to characterize and 
detec t resonances in this context (IContopoulosl . I2OO2I : iLukes-Gerakopoulos et al.l . 



2010 ). In order to compute it we first need to localise the central fixed point uq 



in the Poincare surface of section, see fig. 12.421 Equipped with the set of section 
points of a given trajectory we calculate the angle 9/ = angle(r/, r^+i) between 
the radius vectors of each pair of succeeding section points ut and Rotation 
number v is defined as 

1 ^ 

Xe.. (2.98) 



lim 
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We actually compute the mean of the finite number of values instead of the limit. 
Since the succeeding section points of a given resonant trajectory always skip be- 
tween two islands of stability with a constant orientation, we can infer that the 
overall number of the islands in the chain will be given by the denominator of 
the rotation number when written in a simplest fractional form. Chaotic orbits are 
characterized by stochastic behaviour of both indices, V as well as 

In fig. 12.431 we compute co^/coe and v for a single set of trajectories bound in 
the equatorial potential lobe presented in fig. 12.421 Particles are being launched 
from the equatorial plane differing only in the initial value of radial coordinate 
r(0) separated equidistantly by Ar = 0.0002M, yielding 945 trajectories in to- 
tal. Prominent resonances at V = 1/3 (cOr/coe = 2), v = 2/3 (cOr/coe = 2) and 
V = 1 /2 (co^/coe = 5/2) are apparent and easily identified with the chains of islands 
in fig. 12.421 In the regular regions inhabited by KAM tori the rotation number be- 
haves as non-constant continuous function of r(0) except the intervals of constancy 
corresponding with the resonant chains and two major discontinuities arising from 
our convention of taking 9/ as a positive oriented angle between the radius vectors. 
In the chaotic border regions both indices fluctuate but resonant chains represented 
by constant intervals also appear in this zone. 

Detailed view in fig. 12.441 which zooms the portion of fig. 12.431 reveals that a 
number of faint resonances producing thin Birkhoff chains is present in the chaotic 
zone. Direct detection of these in the surface of section would be rather difficult 
since the readability of the surfaces of section decreases rapidly with increasing 
density of depicted trajectories. In this sense the rotation number appears suitable 
indicator of resonance. 

Resonant chains are in principle detectable in terms of spectral analysis of the 
observed electromagnetic signal. Presence of the Birkhoff chains allows us to dis- 
criminate between perturbed and regular system. Moreover the position and the 
width of the chains reflects other p roperties of the system, strength o f the perturb- 
ing magnetic field for instance. See lLukes-Gerakopoulos et al.l (120 10 ) for a detailed 
discussion of this approach applied to the different type of system of extreme mass 
ratio inspiraling sources emitting gravitational radiation. 

However, the critical assumption of the Poincare-Birkhoff theorem that the 
perturbation of the original integrable system is weak is generally not fulfilled in 
our case. Nevertheless for the trajectories analyzed in figs. l2.42H2.44l the theo- 
rem proved applicable as we managed to identify fragmented tori observed in the 
surface of section with the intervals of constancy of the rotation number at values 
given by simple integer ratios ^, | etc. Therefore we can identify fragmented tori 
in the surfaces of section in fig. l2.42l and in analogous closed equatorial lobes with 
the Birkhoff chains of stability islands anticipated by the theorem. However, we 
did not manage to do so in the case of potential valley in fig. 12.381 neither for the 
off-equatorial lobes, e.g. fig. 12.341 In both cases we observed fragmented tori in 
the section but we could not localise the central fixed point uq which is essential 
for the evaluation of the rotation number. In the potential valleys the central part 
of the section is occupied by the escape corridor and in the off-equatorial lobes the 
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Figure 2.43: Rotation number is compared to the ratio of fundamental frequencies 
Both are computed for the same set of trajectories originating in the equatorial 
plane at different r(0) (equidistantly placed with Ar = 0.0002M). Other parame- 
ters remain fixed: a = 0.5M, E = 1.795, L = 5M, qBo = 2M-\ 9(0) = Gsection = f , 
qQ = 2. Most prominent resonances correspond with v = 1/3 (co^/coe = 2), 
V = 2/3 (co^/coe = 2) and v = 1/2 (co^/coe = 5/2). We note that the resonance 
at r(0) ^ 2.245 M having frequency ratio ^ = 3 is not detected clearly by the ro- 
tation number v which actually oscillates around V = 1/2 instead of being strictly 
constant. Number of the islands comprising the resonant chain is given by the de- 
nominator of the rotation number. Chaotic regions are characterized by a stochastic 
behaviour of the both indicators, V as well as 
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Figure 2.44: Portion of the fig. 12.431 is provided in the high resolution revealing 
the faint resonances producing thin chains of islands whose direct detection in the 
surface of section is almost impossible (see fig. 12.421) . 

structure of the curves in the section is quite different compared to the equatorial 
lobes. Frequency ratios co^/coe corresponding with fragmented tori did not attain 
any simple integer ratio value here. 

We conclude that the Poincare-Birkhoff theorem is only partially applicable 
in our study. Especially the off-equatorial lobes which are actually supported by 
the perturbation (there are no such lobes in the pure Kerr without EM field) fail to 
fulfill the assumption of the perturbation being weak. In the surfaces of section we 
observe the fragmented curves of a different nature than those of resonant chains 
detected in the equatorial lobes. We stress that the issue definitely deserves fur- 
ther attention since the analysis of the fundamental frequencies has observational 
consequences. 



2.10 A magnetic star 



Various types of stars exhibit very strong magnetism. Peculiar main sequence stars 
of Ap and Bp classes may bear l arge-s cale magnetic field of strength ^ 10"^ G 



dBagnulo et all . l2006l : iBorra et all Il982h 



Considerably stronger fields are found 
in the case of degenerate stars. White dwarfs may reach ^ 10^ G dValyavin et al.l . 
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Figure 2.45: For the energy level E = 0.857 we obtain a broad lobe which almost 
touches the surface of the star at r = 4 M. In the upper-left panel, we launch two 
particles with 9(0) = 1.0492, r(0) = 4.75 M. The particle to the left of the star 
starts with w^(0) = and moves chaotically, while the other one with ^'"(O) =0.03 
follows a perfectly regular trajectory. The upper right-panel shows these two types 
of trajectory appear in the surface of section plot. The bottom panels demonstrate 
the difference between the two types in terms of Recurrence Plots. 



20031 ) and magnetars as a magn e tically extremal subclass of neutron stars even 
^ 10^^ G (IDuncan & Thompson! . 1 19921 ). In our context a broad term magnetic 
star encompasses all the above objects. In the following, however, we concentrate 
ourselves mainly on the case of compact stars. 

We describe the gravitational field outside a magnetic star by the Schwarzschild 
metric. 



d^" 



(2.99) 



The associated magnetic field is modeled as a dipole rotating at angular velocity 
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where 



11 =2M^ + 2Mr + r^\og[ 



(2.100) 
(2.101) 

(2.102) 



The related dipole moment M is given by dBakala et al. . 201Clh 



AMhl'\n -2MYl^ Bo 



6M{r^ - M) + 3r^ {r^ - 2M) In (l - 2Mr^^) ' 



(2.103) 



where Bn is the magnetic field at the neutron star equator, is the radius of the star 
surface HI 

We assume eqs. (I2.99M2.100I) to hold outside the star surface (r > n) and 
inside the light cylinder (m'"m^ < 0). We set = 4M as the inner radial boundary 
of the particle motion. As for the light cylinder, the mentioned condition results 
in a relation r^sinQ^Q? = 1 — which implicitly specifies the outer boundary. 
The vector potential (12.100b is valid inside the rigidly corotating magnetospheric 
plasma, which we consider to be an excellent conductor, so that the force-free 
condition F^u^ = holds for the plasma for which = {u',0,0,u'^) and ^ = D.. 

A general formula for the effective potential eq. (12.731) simplifies to the form 



Veff 



3q!M 1lQ.sm^Q 



+ 



8M3 



(2.104) 



1 + 



rsin9 



+ 



3q^M HI sine 
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2.10.1 Motion inside the potential lobes 

Our previous analysis ( Kovaf et al. . 20 id ) revealed a number of distinct types of 
possible topological structures of the effective potential. However it appears that 
the system is not as rich in its dynamical properties. The test particle trajecto- 
ries share some similar features across different classes of the effective potential. 
Therefore, we only present s urveys of particle d ynamics in three exemplary types: 
la, Ila and IIIc (fig. 12.391 see Kovaf et al. ( 2010l ) for the complete review). 



^The existence of extremely compac t stars with ^ 3M is unlikely, but not excluded 
(iBahcall et all 1 19891: IStuchlik et all l2009l) . Most of the realistic equations of state imply a lower 
limit ^ 3 . 5M /Glendenning", ' 1 997^ 



On the other hand, the models of Q- stars do allow a lower 



limit of ^ 2.8M (Miller et al.. 19981 
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Class la lobes grow with energy increasing. Once the level of the equatorial 
saddle point is reached, the lobes merge with each other across the equatorial plane. 
The single merged lobe eventually intersects the surface of the star if the energy 
level is increased sufficiently, letting the particle fall onto the surface. Lobes of Ila 
type also merge via the equatorial plane but in contrast to the first type the merged 
lobe opens toward the light cylinder (beyond which the model becomes invalid). 
Lobes of the class IIIc first open via the off-equatorial saddle points, allowing the 
particles fall onto the star, before the lobes merge through the equatorial plane. 

Now we study the three selected types of the effective potential topology in 
more detail. We are primarily interested whether and how the dynamic regime 
changes across the given range of specific energy E. Especially, we shall address 
what happens with the dynamics when the particle acquires enough energy to cross 
the saddle point. 

The first survey (type la) begins at energy level E = 0.8482, corresponding to 
the closed lobe. In fig. 12.401 we observe that the motion inside the lobe is stable. 
No chaotic properties are detected - neither in Poincare surfaces of section nor in 
the Recurrence Plots. 




Figure 2.46: Regular motion in an off-equatorial lobe of the third type. Parameters 
used: E = 0.99579, L = 6.25382 M, qS^ = 45.87368 L = 6.25382 MM = 
0.011485 M-K Particles are launched from latitude 9(0) = Osection = 1.0492. 

As the energy increases to £ = 0.8485, the symmetrical lobes merge via the 
equatorial plane. By inspecting a number of trajectories in this case we find that 
chaos starts appearing at this point - those particles which notice the gate through 
the equatorial plane always fill the entire allowed region and they move chaotically. 
Nevertheless, there are still such particles which move regularly in one of the two 
parts of merged lobe and never cross the equatorial plane. We can also find transient 
trajectories corresponding to regular motion lasting for some period of time in one 
part of the lobe, followed by chaotic motion over the entire lobe once the particle 
finds and encounters the passage across the equatorial plane. All of the mentioned 
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Figure 2.47: For the energy level E = 0.9962 (other parameters as in fig. 12.461) we 
obtain a large lobe that almost touches the light cylinder and allows the particle to 
fall onto the surface of the star via a narrow passage above and below the equatorial 
plane. The upper left panel shows two trajectories launched from 9(0) = 1.0497, 
w^(0) = 0. One of the particles (starting from r(0) = 61.5 M) follows an unstable 
path and eventually falls on the star surface. On the contrary, the other particle 
(starting from r(0) = 72.5 M) moves regularly and never escapes any given part of 
the lobe. The upper-right panel shows these two kinds of trajectory depicted in the 
Poincare surface of section. Chaotically dispersed points belong to the escaping 
trajectory. The bottom-left panel shows the Recurrence Plot of the escaping par- 
ticle; this plot does not exhibit typical chaotic behavior, although the large-scale 
structures are present. The bottom-right panel presents the Recurrence Plot of sta- 
ble motion. 
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cases are illustrated in fig. 12.411 

Increasing the energy further to £ = 0.857, we obtain a broad potential lobe 
which almost touches the star surface. The situation changes from the previous 
case where the gate connecting the off-equatorial lobes was narrow. Now we do 
not find trajectories which occupy only one part of the lobe, above (or below) the 
equatorial plane. Chaotic trajectories densely filling the entire lobe are typical 
for this setup. We also encounter perfectly regular trajectories forming ribbon- 
like structures spanned between northern and southern borders of the lobe. In 
Poincare sections these appear as regular islands surrounded by a chaotic ocean 
(upper panels of fig. 12.451) . We notice that the RP of this regular trajectory is 
extraordinarily simple and consists of almost perfect diagonal lines (bottom panels 
of fig. 12.451) . Thus its dynamic properties are close to those of a periodic system, 
which is in contrast with the neighboring fully chaotic orbits. 

The second type (class Ila) of the effective potential topology of off-equatorial 
lobes differs from the first one significantly as the lobes do not open towards the 
star when the energy is raised sufficiently. On the contrary, in this case we observe 
that the lobe's boundary touches the light cylinder first. 

The motion in the off-equatorial lobes proves to be regular while the merging 
lobes bring chaos into play. Chaos becomes dominant for broader lobes, however, 
stable regular orbits also persist. The results are similar to those of the first type of 
potential topology discussed above. 

The last analyzed topology of the lobes (class IIIc) differs profoundly from 
the preceding two cases, as can be seen in fig. 12.391 We find that stable motion 
dominates in this setup. This can be verified by comparison with fig. 12.461 

As we further increase the energy level, we obtain more complicated shapes 
of the equipotentials that allow the particle to fall on the star surface. On the 
other hand, opening the outflow gate energetically precedes the merging point of 
both off-equatorial lobes. In fig. I2.47[ we discuss the motion governed by the 
largest possible lobe which almost touches the light cylinder. We observe that 
stable regular orbits are still possible for those particles that do not hit the passage. 

From the above-given discussion we conclude that the motion of charged test 
particles in the off-equatorial lobes allowed by the test field of the rotating magnetic 
dipole on the Schwarzschild background is largely regular. Once the off-equatorial 
lobes merge with each other, chaos may appear. Increasing the energy, the chaotic 
motion becomes typical but, quite surprisingly, very stable orbits also exist under 
these circumstances. 



Chapter 3 

Conclusions 



3.1 Structure of the electro-magnetic field 

In this work we went through various issues concerning the structure of the electro- 
magnetic field which arises from the interplay between the frame-dragging effect 
and the uniform magnetic field with general orientation with respect to the rotation 
axis of the Kerr source. We further generalized the model by allowing the black 
hole to move translationally in a general direction with respect to the magnetic 
background. Components of the electromagnetic tensor F^v describing resulting 
field were given explicitly (in a symbolic way regarding their length) in the terms 
of the former nondrifting solution. Special attention was paid to the comparison of 
various definitions of the electric/magnetic field. We also reviewed the construction 
of three distinct frames attached to the physical observers. 

Before exploring the rich structures arising from the drift and oblique back- 
ground field we first revisited the issue of the expulsion of the aligned magnetic 
field out of the horizon of the extremal Kerr black hole (Meissner effect). Since 
the effect itself has already been discussed thoroughly in the literature we did con- 
centrate on the observer aspect of the problem instead. By comparing alternative 
definitions of the magnetic vector field given in sec. 12.21 combined with the choice 
of the four- velocity profiles presented in sec. 12.41 we came to the conclusion that 
(i) the Meissner effect is observer dependent and (ii) some definitions of the field 
lines do not fit well into the Boyer-Lindquist coordinate system since they artifi- 
cially amplify the effect of the coordinate singularity at the horizon. 

Namely we observed that in the ZAMO tetrad the field does not exhibit the 
Meissner effect while in the frame of freely falling observer (FFOFI) the field is 
expelled. On the other hand in the renormalized field components the expulsion is 
observed for both ZAMO as well as for FFOFI test charges. In coordinate compo- 
nents the Meissner effect also appears but we decide not to use them because the 
coordinate basis is not normalized which causes artificial deformation of the field 
lines. On the other hand the properly normalized physical components appear prob- 
lematic since they amplify the effect of the coordinate singularity at the horizon as 
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mentioned above. We found them not convenient for the use in Boyer-Lindquist 
coordinate system (at least in the region close to the horizon). Asymptotically- 
motivated (AMO) components which are observer independent as they reflect the 
F^v components directly were also employed and the resulting magnetic lines of 
force were identified with the section of the surfaces of the constant magnetic flux 
which represent yet another way to display the field. We note that in the FFOFI 
frame both magnetic and electric fields are expelled out of the horizon in the case 
of the extremal spin. 

Upon introducing the perpendicular component we observe that generally (i) 
the magnetic field is not expelled anymore and (ii) both the electric and magnetic 
fields acquire a tightly layered structure in the narrow zone just above to the hori- 
zon. Structure of the field is surprisingly complex in this region, self-similar pat- 
terns are observed regardless the choice of the observer proving that the layering is 
an intrinsic feature of the field rather than a mere observer effect. 

In the case of BH's translational motion through the aligned field we also ob- 
serve the complex layering of the field which we attribute to the transversal com- 
ponent arising from the Lorentz boost. However, for a sufficiently rapid drift we 
observe a new effect emerging: as the layers transform they give rise to the for- 
mation of the neutral points of both electric and magnetic fields (though not at 
the same location!). The field structure surrounding such point is characterised by 
four distinct domains (bundles of the field lines) divided by two separatrices inter- 
secting at the neutral point. Such a topology is known to result from the separator 
reconnection, a process which has been studied i n the framework of resistive mane- 
tohydrodynamics (MHD), see e.g. IPriest & Forbes (2000. ). In our electro-vacuum 
model, however, it arises entirely from the interaction of the strong gravitational 
field of the rotating BH with the background magnetic field, i.e. it is a mere gravit- 
omagnetic effect. Charged matter injected into the magnetic separator site is prone 
to the acceleration by the electric field since its motion is not affected by the van- 
ishing magnetic field and thus the acceleration is very effective. 

From the astrophysical viewpoint we regard the topological changes which 
the drift causes upon the field structure, especially the formation of the neutral 
points, as our main result demonstrating clearly that the strong gravitation of the 
rotating Kerr source may itself entangle the uniform magnetic field in a surprisingly 
complex way. We suggest that the gravity of the rotating black hole could work as 
a trigger for magnetic reconnection. 



3.2 Motion of charged matter 

We studied the regular and chaotic motion of electrically charged particles near 
a magnetized rotating black hole or a compact star. We employed the method 
of recurrence analysis in the phase space, which allowed us to characterize the 
chaoticness of the system in a quantitative manner. Unlike the method of Poincare 
surfaces, the Recurrence Plots have not yet been widely used to study the chaotic 
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systems in the regime of strong gravity. 

The main motivation for these investigations is the question of whether the mat- 
ter around magnetized compact objects can exhibit chaotic motion, or if instead the 
system is typically regular. One of the main applications of our considerations con- 
cerns the putative envelopes of charged particles enshrouding the central body in a 
form of a fall-back corona, or plasma coronae extending above the accretion disk. 
While we concentrated on the specifications of the RP method in circumstances 
of a relativistic system, the assumed model cannot be considered as any kind of a 
realistic scheme for a genuine corona. We simply imposed a large-scale ordered 
magnetic field acting on particles in a combination with strong gravity. 

Various aspects of charged particle motion were addressed throughout the chap- 
ter. First of all, we investigated the motion in off-equatorial lobes above the horizon 
of a rotating black hole (modeled by Kerr metric equipped with the Wald test field), 
as well as above the surface of a magnetic star (modeled by the Schwarzschild 
metric with the rotating dipolar magnetic field). In both cases we conclude that 
the motion of test particles is regular, which was confirmed for a representative 
number of orbits across the wide range of parameters over all topological types of 
off-equatorial potential wells. This result is somewhat unexpected because the off- 
equatorial orbits require a perturbation to be strong enough (in terms of strength 
of the electromagnetic field), so that it can balance the vertical component of the 
gravitational force. 

Further, we investigated the response of the particle dynamics when the en- 
ergy level E was raised gradually from the potential minimum to values allowing 
cross-equatorial motion. We examined various topological classes of the effective 
potential and came to the conclusion that the cross-equatorial orbits are typically 
chaotic, although very stable regular or bits may also pers i st for a certain intermedi- 
ate energy range. The classical work of iHenon & HeilesI (1 19641) should be recalled 
in this context since it also identifies the energy as a trigger for chaotic motion in 
the analysed simple system. More recently the Heno n-Heiles system was revisited 
in the relativistic context by lVieira & Leteliei (Il996l) . 

We also addressed the question of spin dependence of the stability of motion 
for Kerr black hole in the Wald field. We noticed that this is a rather subtle prob- 
lem. The effective potential is by itself sensitive to the spin value a - hence, we 
had to link the potential value roughly linearly with the energy E to maintain the 
potential lobe at a given position. In other words, we did not find any clear and 
unique indication of the spin dependence of the motion chaoticness. Most tra- 
jectories exhibited regular behavior, which is also in agreement with the previous 
results indicating that motion in off-equatorial lobes is generally regular. On the 
other hand, in the case of the cross-equatorial motion we observed that, for higher 
spins, more chaotic features come into play when compared with the case of slow 
rotation. This trend might be also attributed to simultaneous adjustments of E. 
In other words, it appears impossible to give an unambiguous conclusion about the 
spin dependence of the particles dynamics. Instead, one has to deal with a complex, 
interrelated dependence. 
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In the case of a Kerr black hole immersed in a large-scale magnetic field, we ob- 
served the effect of confinement of particles regularly oscillating around the equa- 
torial plane. Escape of particles from the plane is allowed for a given range of 
initial conditions since the equipotentials do not close; they form an endless axial 
"valley" instead. The escaping trajectories create a narrow, coUimated structure 
parallel to the axis. 



Chapter 4 

Future prospects 



In the following we present our to-do list comprising of the issues which naturally 
arose during the previous study of the electromagnetic fields and charged particle 
dynamics. Most importantly we want to combine ideas of sec. 12.31 and sec. 12.61 In 
other words we plan to investigate ionized particle motion governed by the gen- 
eralized oblique and drifting EM field. Besides that we shall go through several 
rather technical issues related closely to the topic. 

• More general model of gaseous corona. 

We shall enhance our former axisymmetric model by considering oblique 
(misaligned with the rotation axis) magnetic fields in which the central body 
may be uniformly drifting in a general direction. Structure of the electro- 
magnetic field is profoundly enriched and we suppose that similarly the dy- 
namics of the particles will become considerably more complex. We plan 
to discuss the impact of new parameters upon the off-equatorial stable orbits 
and investigate how do they affect the dynamic regime of motion. We will 
try to identify a possible trigger of chaotic dynamics among new parameters. 
Besides standard methods the recurrence analysis will be employed since it 
proved useful in our previous work. 

• Observational consequences. 

We plan to elaborate ideas introduced in sec. l2.9.4l concerning the frequency 
analysis of the off-equatorial orbits. We have seen that fragmented curves we 
observed in the Poincare surfaces of section corresponding with the trajecto- 
ries bound in the closed equatorial lobes may be identified with the Birkhoff 
chains of islands of stability. Such resonant chain is characterized by a sin- 
gle value of a rotation number which is in principle detectable in terms of 
spectral analysis of the observed signal. Presence of the Birkhoff chains al- 
lows us to discriminate between perturbed and regular system. Moreover the 
position and the width of the chains reflects other properties of the system. 
Detailed discus sion of this approach applied to th e different type of system 



may be found in lLukes-Gerakopoulos et al.l (l2010h . However, in our analysis 
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of the off-equatorial trajectories we observed more complicated structures in 
the surfaces of section which do not allow for the straightforward evaluation 
of the rotation number, nor the ratio of fundamental frequencies. Therefore 
we intend to adjust the method for the application to our scenario and infer 
the possible observational consequences for the system of gaseous corona 
we studied theoretically in chap. [2l 

• Magnetic shift of the ISCO. 

Position of the inner edge of the accretion disk is usually identified with the 
marginally stable geodesic orbit r^s (also referred to as innermost stable cir- 
cular orbit, ISCO) whose position is uniq uely determined by the value of 



the black hole spin a dBardeen et al.Lll972h . Common black hole spin mea 



surement met hods are based on t his re lation as they actually determine 



to evaluate a (iMcClintock et al.l . 1201 Ih . In this context we raise the ques 



tion whether the presence of the magnetic field may change th e position of 



ISCO noticeably. Recently a similar problem was addressed by Eakala et al. 



(l201Qb for the case of Schwarzschild source endowed with the dipole mag- 



netic field. An introductory account of the influence of the uniform mag- 
netic field aligned with the symmetry axis of Kerr black hole was brought by 
Prasannal ( 1978 ). We shall discuss the effect of the oblique uniform magnetic 



field around Kerr source upon the marginally stable orbit in detail. 

Application of a new method for the computation of Lyapunov spectra. 
Lyapunov characteristic exponents (LCEs) are the basic indicators of chaos 
which capture the divergent features of the chaotic orbits straightforwardly. 
The classical non-covariant definition of LCEs, however, meets serious dif- 



ficulties in curved spacetimes. Recently IStachowiak & Szydlowskil (120101) 
proposed novel geometrical approach to the computation of the Lyapunov 
spectra which completely avoids the conventional method of solving the 
variational equations to obtain the Lyapunov vectors which are periodically 
Gram-Schmidt orthonormalized along the flow. New algorithm is covari- 
antly formulated and thus seems to be highly convenient for the application 
in general relativistic systems. We plan to implement this method when in- 
specting the dynamics of charged particles. This might be beneficial for both 
the results itselves and also to prove the new method fruitful. 



Appendix A 

Geometrized units 



We use geometrized units instead of SI throughout this work. We set the speed of 
light c, the gravitational constant G, the Boltzmann constant k and the Coulomb 
constant kc = 4^ equal one. 



constant 


SI value 


SI dimension 


geometrized units 


c 


2.998 X 10** 


ms-i 


1 


G 


6.67 X 10-11 


m^s-^kg-i 


1 


k 


1.38 X 10-23 


JK-i 


1 


kc 


8.988 X 10^ 


kgm^s-^C-^ 


1 



Table A.l: Redefined constants in SI and in geometrized units. 



We aim to express arbitrary quantity in the terms of meters (thus it becomes 
"geometrized"). To manage that we construct conversion factor / consisting of 
constants c, G, k and kc whose dimension multiplied by the dimension of the quan- 
tity being converted gives just meters (typically to the power of 1, 2 or —1). To 
make it more clear we convert mass as an example: 

[M]sj = kg while [M]g^^^ = m 
factor needed: [fu] = mkg~^ 
unambiguously: /m = ^ = 7.43 x 10~^^ mkg~^ 
Mgeom = /mMsi = 7.43 X 10-28 Msi mkg-i 
e.g. solar mass: {Mq)^^^^ = 7.43 x IQ-^^ mkg"^ • 1.989 x 10^^ kg = 1472m 
When converting from the geometrized units back to SI we just need to divide 
by the same factor. Factors for basic quantities are given in table A.2. 

Finally we show how to convert specific charge Q which we use in the main 
text often: 



leaving (Q) dimension-less. 

^ V^/geom 
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quantity 


factor 


numerical value 


time 


fi = C 


3.00X lO^ms-i 


mass 


f — G 
JM — 


7.43 X 10-2^ mkg-^ 


charge 


fc = ^ 


8.62X 10-^^mC-^ 


momentum 


f - G 


2.48 X 10"^^ mkg-^m-^s 


angular momentum 


f 

Jm — ^3 


2.48 X 10-^^ m^kg"^m-^s 


energy 


/E — 74 
JT — ^4 


8.26 X 10-45 mJ-i 


temperature 


1.14X 10-^^mK-i 


magnetic induction 


f — 1 

cy kc 


2.87 X 10-1^ m-^T-i 



Table A.2: Conversion factors for basic quantities. 



For example specific charge of the proton is in the geometrized units ex- 
pressed as follows: 

fe)geom = 1-16 X 10^^ • (^p)3j = 1.16 X 10^^ • 9.58 X 10^ = 1.11 X 10^^ 

When dealing with compact objects we often scale all the geometrized quanti- 
ties by the mass of the central body in order to simplify our equations. The mass of 
the object only needs to be specified at the very end of calculations when we need 
to recover actual value of a quantity from it's dimension-less scaled version. 

Scaling of an arbitrary quantity X by the mass of the gravitational source M 
can be formally expressed as follows: 



J geom ^ ^ >' geom, scaled (^M)^ 



geom 



geom 



For instance one reads the particle's orbital proper period from the output of the 
dimension-less equations of motion to be T = 200. Setting the mass of the central 
object as M = 3 Mq we obtain following SI value of the proper period 

(r)geom,scaled(M)geom 200 ■ 3 • 1472 .QzLvin-3. 

(r)si = J = 3X10B =2-94x10 s. 

Angular momentum of the rotating object 5 = a M is also commonly scaled 
by the mass of the object. In the case of the Kerr black hole the spin parameter a 
is restricted to |(^) geom, scaled I ^ 1- ^or example if we consider (highly idealized) 
Sun as a homogeneous sphere of radius = 6.96 x 10^ m rotating with period 
of 25 days (Tq = 2.16 x 10^ s) we arrive at the following value of the scaled spin 
parameter a 

. . (^)geom _ /m(^)si _ fm ( A% MqRI \ _ 
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Kerr solution for the source with given mass and angular momentum would thus 
describe a naked singularity rather than black hole. 

Another example is the product qB of specific charge of test particle (dimension- 
less in geometrized units) and the asymptotic strength of the magnetic field (of 
geometrized dimension m~^) which acts as one of the parameters determining par- 
ticle's trajectory in chap. 12.61 For example setting (^5) geom, scaled = 10 in the 
dimension-less equations would correspond with the following strength of mag- 
netic field B for a given supermassive black hole of mass M = 1 x 10^ Mq if we 
further specify the test particle to be an electron 



geom, scaled seom, scaled 

(qB) 

geom, scaled /m 

(Bjsi - - - 



/b (■W)geom /b (■W)geom(^e)geom, scaled (■W)geom(^e)si /b/c 

10-3 X 10^ 



1 X 106 • 1472- 1.76 X IQii 



= 1.16 X 10-"T= 1.16X IQ-^G. 



Appendix B 

Choice of the integrator 



In this section we shall give some details about a rather technical issue concerning 
the proper choice of the integration scheme which would fit best to our problem. 
In particular we will compare the performance of the symplectic integrator with 
several non- symplectic routines and discuss under which circumstances we should 
choose the symplectic one and when we should switch to some other scheme. We 
will be basically concerned with two crucial aspects - accuracy of the integration 
and CPU time consumption. The latter is generally less critical in our application 
since we are not facing that computationally intensive problem. 

We are dealing with autonomous Hamiltonian system^ whose equations of mo- 
tion form a specific subclass of first order ordinary differential equations (ODEs). 
Two fundamental characteristics of the Hamiltonian flow should be highlighted 

• conservation of the net energy (Hamiltonian) of the system 

• conservation of the symplectic structure CO = dTr^ A dx^. 

In the classical mechanics the natural choice of the generalized coordinates leads 
to the Hamiltonian which may be interpreted as a net energy of the system. This 
is true even for the system of a charged particle in the ex ternal EM field where 



the generalized momenta-dependent potential is introduced (iGoldstein et al.L 12000 , 
chap. 8). Time-independance of the Hamiltonian is thus equivalent to the conser- 
vation of the net energy of the system. In the general rela tivistic version of th is sys- 



tem, however, we employ super-hamiltonian formalism ( iMisner et al.l . ll973L chap. 



21) in which the energy of the particle as a negatively taken time component 
of the canonical momentum E = —71^, is conserved by virtue of the Hamilton's 
equations itselves providing that the super-hamiltonian doesn't depend on the co- 
ordinate time t. On the other hand the value of the super-hamiltonian ^ = ^PjnP^ 



Equations of motion may be equivalently expressed in terms of Lorentz force jMisner et all 

[1921 p. 898) which leads to the set of four second order ODEs. Numerical experiments, however, 
led us to the conclusion that this formulation is computationally less effective compared to the Hamil- 
tonian formalism. Generally for a given numerical scheme with the same parameters (resulting in 
similar accuracy of integration) the integration of Hamilton's equations was roughly two times faster. 
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is by construction equal to —\m^ where m is the rest mass of the particle. Conser- 
vation of the super-hamiltonian in the system is thus equivalent to the conservation 
of the rest mass of the particle. 

By conservation of the symplectic 2-form CO we mean that its components COo^p 
in the basis (dr(?i),dr(?i),d9(?i),d(p(?i)),d7r^(?i),d7i^(?i),d7i0(?i),d7i(p(?i)) do not 
change during the evolution of the system and for arbitrary value of the affine 
parameter X (i.e. at each point of the phase space trajectory) we obtain 



-I 




(B.l) 



where I stands for the four-dimensional identity submatrix and is null subma- 
trix of the same dimension. Conservation of the symplectic structure expresses in 
the abstract geometrical language the fact that th e evolution of the system is gov- 
erned by the Hamilton's canonical equations. See Arnold ( 1989 ) for details on the 



geometric formulation of the Hamiltonian dynamics. 

It would be highly desirable to use such integration scheme which would con- 
serve both quantities which are conserved by the original system. It appears, how- 
ever, that this is not possible for non-integrable systems and one has to decide 
whether he employs the scheme which conserves energy or rather the integrator 
which keeps symplectic structure. The latter are referred to as symplectic inte- 
grators and by many accounts provi de most reliable results in numerical studies 



involving Hamiltonian systems. See lYoshidal (Il993h for a comprehensive review 
on symplectic methods. 

We list all the schemes we employ in this survey specifying their basic prop- 
erties. We shall actually compare one symplectic method with several standard 
integrators. Code names we use for the schemes are those which denote the rou- 
tines in the MATLAB system. 

• GLS - Gauss-Legendre symplectic solver, ^-stage implicit Runge-Kutta (RK) 
method, crucial control parameter: stepsize h 

• ODE87 - Dormand-Prince 8th - 7th order explicit RK scheme, the most 
precise RK method (local error of order 0{h^)), adaptive stepsize - RelTol 
is set to control local truncation error 

• ODEl 13 - multistep Adams-Bashforth-Moulton solver, based on the predic- 
tor-corrector method (PECE), RelTol is set 

• ODE45 - Dormand-Prince seven stage 5th-4th order method of explicit RK 
family, adaptive stepsize, default integration method in MATLAB and GNU 
OCTAVE, error is controlled by RelTol 

Apart from ODEl 13 all other routines are single-step (Runge-Kutta like) methods 
which means that they express the value of the solution in the next step in terms 
of a single preceding step. They may be related explicitely or implicitly. Multistep 
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methods in contrast employ more preceding steps to calculate the solution at the 
succeeding point. RelTol is a parameter which specifies the highest allowed relative 
error in each step of integration (local truncation error) when the adaptive stepsize 
methods are used. In the case of exceeding the RelTol the stepsize is reduced 
automatically to decrease the error. 




Figure B.l: Regular trajectory of charged test particle (qQ = 1, L = 6 M and E = 
1.6) on the Kerr background (<2 = 0.9 M) with Wald magnetic field (qBo = IM"^). 
Particle is launched at r(0) = 3.68, 9(0) = 1.18M with u'{0) = 0. 



We comment that for general non- separable Hamiltonians only implicit sym- 
plectic schemes may be found. Explicit methods e xist for sepa rable Hamiltonians 
and for some special forms of non-separable ones (IChinl . l2009h . Besides other im- 
plications of the usage of the implicit methods we note that they necessarily involve 
some type of iterative scheme which is typically of a Newton's type and thus re- 
quires to supply Jacobian of the right hand sides of the equations of motion which 
is the Hessian matrix of the second derivatives of the super-hamiltonian !H in our 
case. 

Another inconvenience connected with the symplectic methods is their fail- 
ure to co nserve the symplectic structure once the adaptive stepsize method would 
be used (iSkeel & Geai . Il992h . Therefore the stepsize has to be set rigidly for a 
given integration segment when using symplectic method. Several workarounds 
have been suggested to combine benefits of symplectic so lvers and vari able step- 
size algorithms - e.g. Hairer's symplectic meta-algorithm (lHaireij . ll997n which is, 
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affine parameter X 

Figure B.2: Comparison of the integrators in the case of regular trajectory. Sym- 
plectic GLS provides the most reliable results for X> 10^. Bottom panel shows 
that besides secular drift in energy (artificial excitation or dumping of the system; 
plot shows absolute values, however) it also oscillates on the short time scale. 
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however, only applicable to the separable Hamiltonians. In our context one would 
considerably suffer from the fixed timestep only in the case of highly eccentric 
orbits. 



integrator 


A\E\/\E\ 


^comp [h] 


RelTol 


stepsize h 


GLS 


« 10-1'' 


14 


N/A 


0.25 


ODE87 


«10-9 


14 


10-14 


adaptive 


ODE113 


« 10-3 


1/3 


10-14 


adaptive 


ODE113 


« 10-3 


1/4 


10-6 


adaptive 


ODE45 


W 10-3 


1/4 


10-14 


adaptive 



Table B.l: Comparison of the performance of several integration schemes for the 
regular trajectory integrated up to ?i = 4 x 10^ (see fig. IB. 21) . 




Figure B.3: Chaotic trajectory of charged test particle (qQ = 1, L = 6 M and E = 
1.8) on the Kerr background (a = 0.9 M) with Wald magnetic field (qBo = 1M~^). 
Particle is launched at r(0) = 3.68 M, 9(0) = 1.18 with w"(0) = 0. 

First we integrate the cross-equatorial regular trajectory fig. IB. 11 Comparison 
of the performance of the integrators is plotted in fig. IB. 21 We plot relative devia- 
tion of the particle's specific energy E from its initial value rather than the error in 
super-hamiltonian because the discussion of motion in chap. 12.61 was mostly held 
in terms of E whose impact upon the trajectory is thus more familiar to us. We 



APPENDIX B. CHOICE OF THE INTEGRATOR 



101 




0.5 1 1.5 2 2.5 3 3.5 4 

a*fine parameter ^„ ^ 




1000 2000 3000 4000 5000 6000 7000 8000 

affine parameter X 



Figure B.4: Comparison of the integrators in the case of chaotic trajectory. For 
'k>5x 10^ the GLS dominates in accuracy over other schemes with the difference 
rising steadily. In the upper panel we compare ODEllS's outcome for two distinct 
values of RelTol parameter. ODE45 is not shown to avoid overlapping of its plot 
with ODE133 curves. 
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calculate the current value of E from the super-hamiltonian !H , while the value of 
7Zt remains truly constant regardless the integrator since the Hamilton's equation 
for its evolution is simply diZt/dX = 0. 

Stepsize of GLS is set in such a way that the integration consumes roughly 
the same amount of the CPU time as it does for ODE87 with RelTol = 10"^"^ 
to make the results comparable. The global accuracy of the GLS solver could 
be further increased by reducing the stepsize while decreasing the RelTol hardly 
improves the secular accuracy of non-symplectic methods here (we have compared 
RelTol = 10-^ and RelTol = 10-^"^ results for ODE113 obtaining global errors of 
the same orders in both cases). 

We observe that the error of GLS rises steeply at the beginning and ODE87 
is considerably better for some amount of time. However then the error of GLS 
almost saturates while ODE87's error keeps growing significantly. For ?i > 10^ 
which corresponds to ^ 1000 revolutions around the center^ the GLS scheme 
becomes more accurate than ODE87 with the difference further rising steadily. We 
conclude that in the case of regular trajectory ODE87 is appropriate for short-term 
accurate integration and GLS for any longer accurate integrations. On the other 
hand for fast, though inaccurate computations one employs ODE113 on all time 
scales. 



integrator 


\AE\/ 


\E\ 


^comp [h] 


RelTol 


stepsize h 


GLS 


lo- 


-9 


14 


N/A 


0.25 


ODE87 


ss 10" 


-6 


14 


10-14 


adaptive 


ODE113 


« lo- 


-3 


1/6 


10-14 


adaptive 


ODE113 


ss 10" 


-3 


1/6 


10"^ 


adaptive 


ODE45 


10" 


-3 


1/2 


10-14 


adaptive 



Table B.2: Comparison of the performance of several integration schemes for the 
chaotic trajectory integrated up to ?i = 4 x 10^ (see fig. lB.4l) . 

In the case of the chaotic trajectory fig. IB.4I the dynamics changes in favor 
of symplectic solver GLS. In fig. IB. 31 we observe that in this case the symplectic 
scheme is superior to the others in even more convincing manner than it was in 
the regular case. Although the initial phase when the error induced by GLS rises 
more steeply than that of ODE87 is also present, it turns over very quickly and for 
?i > 5 X 10^ 50 azimuthal revolutions) the GLS turns out to be more accurate. 
The difference then rises much faster compared to the regular case. 

Experiments with ODE113 show that here we observe distinct (thought not 
sharply) errors by changing the RelTol. Difference of eight orders of magnitude 
in RelTol resulted in roughly one order difference in global error. We also note 
that chaotic regime induces disorder in short-time oscillations of the global error 

^For instance for M = lO^M© the azimuthal proper period of a given particle reads T(p^ lO^s in 

SI. 
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(see bottom panel of fig. IB.4I) . We summerize that the chaotic regime accents the 
supremacy of the symplectic scheme which is to be applied on all time scale here 
(except very short integrations where ODE87 dominates) to obtain the most accu- 
rate results. For fast though inaccurate calculation one would switch to ODE113 
as before. 



GLS, 1^=3 Jt 10^ 0DEe7, K=3y^^(f 




s.Gaaa 2^01 2^012 2^014 2.Gciia 2.eais 2^02 2^022 2.Gaaa 2^01 2.^012 2.^qi4 z«iq 2.S113 2.502 2.^022 

r[M] r(M| 



Figure B.5: We show how the accuracy of the integration crucially affects the 
appearance of the Poincare surfaces of section of a single regular trajectory with 
qQ = 1.76, L = 4.02 M and ^ = 1.619855 on the Kerr background a = 0.55 M 
with Wald magnetic field qBo = 1.92 M'K Particle is launched at r(0) = 2.5012, 
9(0) = 1.0447M with u'{0) = 0. We distinguish u^>0 (black point) from <0 
(red point) in the surfaces of section. 

From a practical point of view we demand high accuracy of the long-term inte- 
gration when constructing Poincare surfaces of section. By theory the intersection 
points with regular trajectory form one-dimensional curve in the section plane. In 
fig. IB. 51 we observe, however, that the points may be dispersed over the consid- 
erable area if the global error in energy rises causing artificial excitation/dumping 
of the system. Symplectic GLS provides the most reliable outcome, with ODE87 
the curve is blurred significantly but the interpretation remains unambiguous. With 
ODEl 13 the curve is further blurred and using ODE45 solver we obtain completely 
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unreliable outcome which could easily lead to incorrect interpretation of a trajec- 
tory as a chaotic one. We note that we intentionally chose such trajectory which 
is highly sensitive to the relative errors in dynamic quantities since it itself spans 
small range of coordinate and momenta values. 

We confirm that the symplectic integrators are the method of choice in the 
case of long-term integration of the Hamiltonian system which in our case con- 
sists of a charged test particle orbiting around the Kerr black hole with stationary 
and axisymmetric electromagnetic test field. Its supremacy over non-symplectic 
methods is even more apparent in the case of chaotic orbits, where the global accu- 
racy of non-symplectic methods decreases rapidly. The accuracy of the symplectic 
integrator could be further increased by reducing the stepsize (at the cost of the 
computational time). On the other hand the performance of the non-symplectic 
solvers is not considerably affected by changing the RelTol parameter across the 
wide range of the values. Once the integrator doesn't fit the problem (= is not sym- 
plectic) there is no effective way to control the global error and even the extremely 
small local truncation errors do not ensure reliable outcome on a long time scale. 
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Here we briefly introduce the software tool vfexplorerl we programmed in Matlab 
system in order to survey electric and magnetic fields studied in chap. [2l We set up 
simple graphic user interface (GUI) for the sake of effective exploration of complex 
field structures. 

Basic features of vfexplorer2 involve interactive zooming, step- wise navigation 
and track logging allowing the user to return to any previous step of the exploration. 
Moreover, he can export the plot to the graphic file directly at any moment. Alter- 
natively he may redirect the matrix out put describing given portion of the vector 
field to the LIC toolbox dShambd . l2005h which encodes the field structure into the 
texture resembling iron fillings (see e.g. fig. l2.17l) 

In fig. lC.ll we present two screenshots illustrating the survey of magnetic field 
around a drifting Kerr BH (see sec. 12.31 for details). Interactive zooming feature 
proves very useful when searching for a particular topological pattern (e.g. separa- 
tor null points which we studied in sec. 12.5.21) in the complex field structure. 

Currently the vfexplorerl GUI also handles the problem specific parameters 
such as spin of the BH, velocity of the drift etc. Nevertheless, it can be adopted to 
explore any other vector field quite easily. Besides that we plan to include another 
features, especially we intend to encode the field strength into the color scale of the 
background of the field lines. 
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Figure C.l: Survey of the complex magnetic field above the drifting extremal Kerr 
BH using the vfexplorer2 GUI. Tightly layered structure near the horizon observed 
in the upper panel is zoomed in the bottom one. 
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